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SCALAR CONSERVATION LAWS WITH ROUGH FLUX AND STOCHASTIC 

FORCING 

MARTINA HOFMANOVA 


Abstract. In this paper, we study scalar conservation laws where the flux is driven by a geometric 
Holder p-rough path for some p G (2, 3) and the forcing is given by an Ito stochastic integral driven 
by a Brownian motion. In particular, we derive the corresponding kinetic formulation and define an 
appropriate notion of kinetic solution. In this context, we are able to establish well-posedness, i.e. 
existence, uniqueness and the L 1 -contraction property that leads to continuous dependence on initial 
condition. Our approach combines tools from rough path analysis, stochastic analysis and theory of 
kinetic solutions for conservation laws. As an application, this allows to cover the case of flux driven 
for instance by another (independent) Brownian motion enhanced with Levy’s stochastic area. 


1. Introduction 

The goal of the present paper is to develop a well-posedness theory for the following scalar rough 
conservation law 

du + div (A(x, u)) d z = g{x , u) d W, t G (0, T), x G 

u(0) = u 0 , 

where 0 = (z 1 ,.z M ) is a deterministic rough driving signal and W = (W 1 ,..., W K ) is a Wiener 
process and the stochastic integral is understood in the Ito sense. The coefficients A : I w xR —» R NxM , 
g : xl-> R a " satisfy a sufficient regularity assumption introduced in Section 2. The above equation 

can be rewritten by using the Einstein summation convention as follows 

du + d Xi u)) d z-’ = gk(x, u) d W k , t e (0, T), x € R N , 

u( 0) = u 0 - 

As an application of our analysis, one can replace z, for instance, by another Brownian motion B , 
which is independent of W, and give meaning to 

du + div (A(x, u)) o d B = g(x, u ) d W, t £ (0, T), x £ R^, 

(L2) u(0) = uo- 

Conservation laws and related equations have been paid an increasing attention lately and have 
become a very active area of research, counting nowadays quite a number of results for deterministic 
and stochastic setting, that is for conservation laws either of the form 

(1.3) d t u + div (A(u)) = 0 
(see [3], [4], [24], [28], [32], [33], [37], [38]) or 

(1.4) du + div (A(u))dt = g(x, u)dW 


Date : March 24, 2016. 

2010 Mathematics Subject Classification. 60H15, 35R60, 35L65. 

Key words and phrases, scalar conservation laws, rough paths, kinetic formulation, kinetic solution, BGK approxi¬ 
mation, method of characteristics. 


1 



2 


MARTINA HOFMANOVA 


where the Ito stochastic forcing is driven by a finite- or infinite-dimensional Wiener process (see [2], 
[7], [11], [12], [13], [16], [21], [23], [27], [41], [42]). Degenerate parabolic PDEs were studied in [4], [8] 
and in the stochastic setting in [1], [10], [22], 

Since the theory of rough paths viewed as a tool that allows for deterministic treatment of stochastic 
differential equations has been of growing interest recently, several attempts have already been made 
to extend this theory to conservation laws as well. First, Lions, Perthame and Souganidis (see [30], 
[31]) developed a pathwise approach for 

d u + div (A(x, u)) o d W = 0 

where W is a continuous real-valued signal and o stands for the Stratonovich product in the Brownian 
case, then Friz and Gess (see [17]) studied 

du + div f(t, x, u)dt = F(t, x, u)dt + Ak(x, u, S7u)dz k 

where is affine linear in u and Vu and z = (z 1 ,..., z K ) is a rough driving signal and Gess and 
Souganidis [20] considered 

du + div {A(x, w))dz = 0 

where 2 = (z 1 ,..., z M ) is a geometric a-Holder path and in [19] they studied the long-time behavior 
for z being a Brownian motion. 

In order to find a suitable concept of solution for problems of the form (1.1), it was observed already 
some time ago that, on the one hand, classical C 1 solutions do not exist in general and, on the other 
hand, weak or distributional solutions lack uniqueness. The first claim is a consequence of the fact that 
any smooth solution has to be constant along characteristic lines, which can intersect in finite time 
(even in the case of smooth data) and shocks can be produced. The second claim demonstrates the 
inconvenience that often appears in the study of PDEs and SPDEs: the usual way of weakening the 
equation leads to the occurrence of nonphysical solutions and therefore additional assumptions need to 
be imposed in order to select the physically relevant ones and to ensure uniqueness. Hence one needs 
to find some balance that allows to establish existence of a unique (physically reasonable) solution. 

Towards this end, we pursue the ideas of kinetic approach, a concept of solution that was first 
introduced by Lions, Perthame, Tadmor [33] for deterministic hyperbolic conservation laws and further 
studied in [3], [8], [24], [32], [33], [38], [37]. This direction also appears in several works on stochastic 
conservation laws and degenerate parabolic SPDEs, see [10], [11], [12], [13], [21], [22] and in the (rough) 
pathwise works [19], [20], [30], [31]. In comparison to the notion of entropy solution introduced by 
Kruzkov [28] and further developed e.g. in [2], [4], [16], [27], [42], kinetic solutions are more general 
in the sense that they are well defined even in situations when neither the original conservation law 
nor the corresponding entropy inequalities can be understood in the sense of distributions which is 
part of the definition of entropy solution. Usually this happens due to lack of integrability of the flux 
and entropy-flux terms, e.g. A(u) ^ Lj . Therefore, further assumptions on initial data or the flux 
function A are in place in order to overcome this issue and remain in the entropy setting. It will 
be seen later on that no such restrictions are necessary in the kinetic approach as the equation that 
is actually to be solved - the so-called kinetic formulation, see (2.3) - is in fact linear. In addition, 
various proofs simplify as methods for linear PDEs are available. 

Let us now shortly present the main ideas of our approach. Apart from the above mentioned 
difficulties there is another one that originates in the low regularity of driving signals and solution. 
Namely, the corresponding rough integrals are not well defined so we present a formulation, see (2.5), 
that does not include any rough path driven terms and therefore provides a suitable notion of kinetic 
solution in this context. To this end, we adapt the ideas of [30], [31], [20] where the authors introduced 
a method of modified test functions to eliminate the rough terms. Our method then combines this 
approach with the ideas for stochastic conservation laws treated in [12], [21]. As usual for this class 
of problems, we define a second notion of solution - a generalized kinetic solution which, roughly 
speaking, takes values in the set of Young measures. The general idea is that in order to get existence of 
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such a solution, weak convergence (in some L p space) of the corresponding approximations is sufficient 
which allows for an easier proof. The key result is then the Reduction Theorem that states that 
any generalized kinetic solution is actually a kinetic one, that is, the Young measure at hand is a 
parametrized Dirac mass. 

Concerning the existence part, we make use of the so-called Bhatnagar-Gross-Krook approximation 
(BGK for short) which allows to describe the conservation law as the hydrodynamic limit of the 
corresponding BGK model, as the microscopic scale e vanishes. This is nowadays a standard tool 
in the deterministic setting where the literature is quite extensive (see [3], [24], [32], [33], [38], [37]). 
Even though the general concept is analogous in the stochastic case, the techniques are significantly 
different, the result was established recently by the author [21]. To be more precise, the key point is 
to solve the corresponding characteristic system which for the general case of (1.1) reads as follows 

d <Pt = -d Xi A ij (<p t )dz :i +g k (vt)dW k , 

dtp\ = duAij&t) dz 3 , i = l,...,N. 

Already at this level one can understand the difficulties coming from the complex structure of (1.1). 
Namely, in the case of (1.3) the characteristic system reduces to a set of independent equations (note 
that the flux function A is independent of the space variable) 

d ‘fit = 0 

d <p\ = duAtj (ip°)dt, i = l,...,N, 

which can be solved immediately and as a consequence the method simplifies. For the stochastic case 
(1.4) we obtain 

d <fit = g k {vt)dW k 

M 

d At = ^duAij^dt, i = l,...,N, 
j=i 

hence the first coordinate of the characteristic curve is governed by an SDE and further difficulties 
arise due to randomness (see [21]). If A is also x-dependent we observe an additional term in the first 
equation of the characteristic system (1.5), however, let us point out that there is a major difference 
between this and the term coming from the forcing. In particular, if g = 0 then the flow of diffeo- 
morphisms generated by (1.5) is volume preserving as can be seen easily by calculating divergence of 
the corresponding vector field. This does not hold true anymore if forcing in nonconservative form is 
present in the equation, unless g is independent of the solution, i.e. we consider additive noise in (1.1). 

The method of BGK approximation has another advantage especially when dealing with rough 
path driven conservation laws. As already mentioned above, the problem boils down to solving the 
characteristic system (1.5), i.e. an ordinary (stochastic, rough) differential equation and the theory 
for those problems is well-established, unlike the theory for rough path driven equations in infinite 
dimension, i.e. rough PDEs. Furthermore, the BGK model also provides an explicit formula for the 
approximate solutions and therefore the necessary estimates independent of the microscopic scale £ 
come rather naturally. 

The exposition is organized as follows. In Section 2, we introduce the basic setting, define the 
notion of kinetic solution and state our main result, Theorem 2.8. In order to make the paper more 
self-contained, Section 3 provides a brief overview of the relevant concepts from rough path theory. 
Section 4 is devoted to the proof of uniqueness, reduction of a generalized kinetic solution to a kinetic 
one and the ZA-contraction property, Theorem 4.6 and Corollary 4.7. The remainder of the paper 
deals with the existence part of Theorem 2.8 which is divided into several parts and finally established 
through Theorem 5.9. Subsections 5.1 and 5.2 give a rough-pathwise existence of unique solutions to 
the BGK approximation. In these two sections we work with a fixed realization of the driving signals 
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or, to be more precise, a fixed realization of a joint lift of (t,z,W). In Subsection 5.3, the stochastic 
approach is resumed and we pass to the limit and obtain a kinetic solution to ( 1 . 1 ). 


2. Definitions and the main result 


Let us now introduce the precise setting of (1.1). We work on a finite time interval [0 ,T\, T > 0 and 
on the whole space concerning the space variable x £ R w . Throughout the paper, we assume that z 
can be lifted to a geometric Holder p-rough path for some p £ (2,3) and denote its lift by z = (z 1 , z 2 ). 
Next we consider the following joint lift of (z, W): we define A = (A 1 , A 2 ) by 


( 2 . 1 ) 


4, 

if i £ {1,.. 

-,M}, 



wt M , 

if i £ {M 4 

"I,- 

..,M + K}, 





if 

i,j £ {1,..., M}, 


fowr M 

0 d wl~ M 


if 

i,j €{M + 1,...,M + K}, 


Jo4dW? 

-M 


if 

i £ {!) • ■ ■ j £ {M + 1. 


Awt M 

~ fo 4 dW r~ 

-M 

if 

i £ {M + 1,...,M + K},j 

£ {1,..., M}. 


It was shown in [14] that such a stochastic process exists and could be considered as the canonical joint 
lift of z and W. Note that although in the original equation (1.1) we consider Ito stochastic integral, 
the lift of W used in the construction of A above corresponds to the Stratonovich version. 

Regarding the coefficients in (1.1), let us fix the following notation. Let 

a = (ay) = (d u A i:j ) : R N x R —> R NxM , 

b = (bj) = (div x A.j) : R n x R —» R M , 


and assume that a, b £ Lip 7+2 and g £ Lip 7+,i for some 7 > p. Here we adopt the notation of [18, 
Definition 10.2], namely, a mapping V : R e —► R d belongs to Lip^ provided it is bounded, [_/3J -times 
continuously differentiable with bounded derivatives of all orders and its |_/3J th derivative is {/3}-Holder 
continuous. 

Furthermore,we suppose that 

( 2 . 2 ) b{x, 0 ) = 0 g{x, 0 ) = 0 Vx £ 1 ^, 


and denote 


K 

G 2 {x,0 = Y J \9k(.x,0\\ VzeR^eR. 

fc=1 


2.1. Notations. We adopt the following notations. The brackets (■, •) are used to denote the duality 
between the space of distributions over R w x R and C* (R ,v x R). We denote similarly the integral 

(/, h) = I [ f(x,€)h(x,Z)dxd£, f £ L P (R N x M), h £ L q (R N x M), 

J R w J K 

where p, q £ [l,oo] are conjugate exponents. By _Mf,([0,T) x 1^ x R) we denote the space of Borel 
measures over [0,T) x M. N x R and A4^([0,T) x M. N x R) are then nonnegative Borel measures. We 
also use the shorthand 

n{<j>) = f (j>(t,x,£)dn(t, x,£), n £ A4b([0, T) x R w x R), cj> £ C c ([0, T) x R w x R). 

J[0,T)xR Ar xK 

and Co([0, T] xR ff xR) denotes the space of continuous functions [0, T\ x R w xR that vanish at infinity, 
i.e. for large ( x, £). The differential operators gradient V and divergence div are (unless otherwise 
stated) understood with respect to the space variable x. 
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2.2. Definitions. As the next step, let us introduce the kinetic formulation of (1.1) as well as the basic 
definitions concerning the notion of kinetic solution. The motivation behind this approach is given by 
the nonexistence of a strong solution and, on the other hand, the nonuniqueness of weak solutions, 
even in simple cases. The idea is to establish an additional criterion - the kinetic formulation - which 
is automatically satisfied by any strong solution to (1.1) (in case it exists) and which permits to ensure 
the well-posedness. 

We start with the definition of kinetic measure. 


Definition 2.1 (Kinetic measure). A mapping to from fl to ([0, T] xl^xR), the set of nonnegative 
bounded measures over [0,T] x R N x R, is said to be a kinetic measure provided 

(i) m is measurable in the following sense: for each <f> £ (7o([0,T] x x R), the mapping 
m(0) : S2 —>• R is measurable, 

(ii) the mapping 

dm(t, x, £) : Q —>• R 


is measurable and 


(iii) for any <p £ Cq 


'[0,T)xK N xt 


E 


dm(t, x, £) < oo, 


J[0,T)xl w xI 

x R), t eA TO(l[ 0j t]</>) is progressively measurable. 


Formally speaking, the kinetic formulation corresponding to the conservation law at hand is given 
as follows 


(2.3) 


dF + X7F ■ adz — d^F bdz= -d^FgdW + ( G 2 d i F ) d t + <%to, 

^(0) = F 0 , 


where F = l u >£ and to is a kinetic measure 1 . However, since the expected regularity of solutions is 
low and consequently the rough path driven integrals are not well defined, it is necessary to define 
a suitable notion of weak solution to this problem. This leads us to the notion of kinetic solution 
to rough path driven conservation laws that we introduce in this work. Note that it is a consistent 
extension of the corresponding notion of kinetic solution for the case of a smooth driving signal z, for 
further discussion on this subject we refer the reader to Subsection 2.3. 

Definition 2.2 (Kinetic solution). Let uq £ L 1 n L 2 (fi x R^). Then a progressively measurable 

u £ L 2 (H;L 2 (0,T; L 2 (R w ))) 


satisfying 


(2.4) 


Eesssup ||u(t)||ii < C 

o <t<T 


is said to be a kinetic solution to (1.1) with initial datum uq provided there exists a kinetic measure 
to such that the pair (F = l M >f,m) satisfies, for all </> £ x R) and a £ C^([0,T)), P-a.s., 


(2.5) 


[ { F {t)A{6t))dta{t) dt + (F 0 , (f>)a{ 0) 

Jo 

= [ {9^F(t)g, 4>(9 t ))a(t) dW — f (<%(G 2 d ? F(t)), 0(6> t ))a(t) dt + m(ad^(e t )), 
Jo z Jo 


^Here it is a function of (cj,t,x) so F(u,t, x, £) = l u (u,t,x)>£ i s well-defined and regarded as a function of four 
variables x,£). 
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where 9 = ( 9 °, 9 X ) is the inverse flow corresponding to 

d7r° = —b(-Kt) dz 
d7r“ = a(7r t ) dz. 


( 2 . 6 ) 


To be more precise, it follows from [18, Proposition 11.11] that under our assumptions (2.6) possesses 
a unique solution 7r that defines a flow of C 2 -diffeomorphisms. We denote by 7r Sit (x,£) the solution 
of (2.6) starting from (x,f) at time s. To simplify the notation, we write 7r t instead of and we 
denote the corresponding inverse flow by 9. Then 9t }S = is the unique solution to the time-reversed 
problem 

(27) d9 t,s = -b(9t, s ) dz s , 

d 9f s = a(6» tjS ) dz s , 

where z s (-) = z(s — •) is the time-reversed path to z. We point out that the flow n as well as 9 is 
volume preserving as can be seen easily by calculating divergence of the corresponding vector field 
and recalling the fact that divergence free vector fields generate volume preserving flows. Thus the 
Jacobian of 7r Sjt satisfies: J7r S[t = 1 and similarly for 9. 

Note that with a classical argument of separability, the set of full probability where (2.5) holds true 
does not depend on the particular choice of test functions (f> , a. 

We proceed with a reminder of Young measures and the related definition of kinetic function that 
will eventually lead to the notion of generalized kinetic solution, see Definition 2.7. The concept of 
Young measures was developed in [43] as a technical tool for describing composite limits of nonlinear 
functions with weakly convergent sequences, for further reading we refer the reader e.g. to [36]. 

In what follows, we denote by 'Pi(R) the set of probability measures on R. 


Definition 2.3 (Young measure). Let (A, A) be a a-finite measure space. A mapping v : A — » 'Pi(R) 
is called a Young measure provided it is weakly measurable, that is, for all <fi £ C'b(R) the mapping 
z ifrom A to R is measurable. A Young measure v is said to vanish at infinity if 


|£|dzA,(£)dA(2) < oo. 


Definition 2.4 (Kinetic function). Let (A, A) be a u-finite measure space. A measurable function 
/ : X x R —>• [0,1] is called a kinetic function on X if there exists a Young measure ronl that 
vanishes at infinity such that for a.e. z £ X and for all £ £ R 


f(z,0 = i" s (£,oo). 


Lemma 2.5 (Compactness of Young measures). Let (A, A) be a a-finite measure space such that 
L 1 (X) is separable. Let ( v n ) be a sequence of Young measures on X such that for some p £ [1, oo) 

(2.8) sup / [ \f\ p di>z(£)d\(z) < oo. 

n£N JX J R 

Then there exists a Young measure v on X and a subsequence still denoted by ( v n ) such that for all 
h £ L 1 (A) and all <j> £ C'b(R) 

lim [ h(z) [ <K£)di/”(£)dA (z) = [ h{z) [ 0(£) di/ z (£) dA(z) 

n_> °° J x Jr Jx Jr 

Moreover, if f n , n £ N, are the kinetic functions corresponding to v n , n £ N, such that (2.8) holds 
true, then there exists a kinetic function f (which correponds to the Young measure v whose existence 
was ensured by the first part of the statement) and a subsequence still denoted by ( f n ) such that 


L°°(X x R). 
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Proof. Various results of this form are somewhat classical in the literature, a proof for the case of 
(X, A) being a finite measure space can be found in [12, Theorem 5, Corollary 6], however, one can 
actually observe that this additional assumption is not used in the proof and therefore the same proof 
applies to our setting of (X, A) being a cr-finite measure space. □ 


Remark 2.6. If / : X xR -A [0,1] is a kinetic function corresponding to the Young measure v satisfying, 
for some p G (1, oo), 


[ [ |£| p du 2 (£) dA(z) < oo, 


JX J R 

then we denote by Xf the function defined by Xf( z i 0 = f{ z i £) — lo>{- Contrary to /, this modification 
is integrable on E x R whenever E C X with A (E) < oo. Indeed, 


hence 


which implies 


Xf 0,0 


/(-oo,£] C < 0) 

f(£,oa) ^ Vz ' € — 


ICI P / lx/(z,0|dA(z) < [ f \C\ p dv z (Qd\(z) <oo 

JX JXJR 


JeL l X/ ( 2 ’^)l d ^ dA ( Z ) - J R 1 + 1 |g|P d ^ + J x J R ICI P du 2 (C)dA(z)^ < oo. 

Besides, if / is at equilibrium, that is, there exists u G L 1 {X) such that / = l u >^ and v = <5„={ then 
we will rather write Xu instead of Xf and it holds true that 

/ |x«(*)(£)|d£ = |u(z)|. 

Jr 


Definition 2.7 (Generalized kinetic solution). Let Fq : tt x R w x 1 A [0,1] be a kinetic function. A 
progressively measurable function F : XI x [0, T\ x x K A [0,1] is called a generalized solution to 
(1.1) with initial datum Fq if F(t) is a kinetic function for a.e. t G [0,T], there exists C > 0 such that 

(2.9) E css sup [ [ |£|du tiX (£)da; + E f f f |£| 2 di/ t ,x(€) dx < C, 

o <t<T J i» Jr Jo Jrn J r 

where v = —d^F, and if there exists a kinetic measure m such that (2.5) holds true. 


2.3. Kinetic solutions in the case of smooth driver z. For the sake of completeness, let us now 
present a formal derivation of (2.3) in case of sufficiently smooth solution to the conservation law 
(1.1) driven by a smooth path z and a Wiener process W. We denote such a solution by u and note 
that it satisfies (1.1) pointwise in x G 1^. Let us fix x G R w and derive the equation satisfied by 
lu(-,x)>£ which is understood as a distribution in £. Towards this end, we denote by (•, -)j the duality 
between the space of distributions over R and C£°(R) and observe that if </>(£) = f^_ <j> i(C)d£ for 

some <j>i G C£°(R), then 

(1 = </>(«) - </>(-oo) = 

Therefore using the Ito formula, we obtain 

d {1„>£, = —(j>\(u) div (A(x, u))dz + cj>i(u)g(x, u)dW + -cj)' 1 {u)G 2 {x, u)dt 

= —(j>i(u)a(x,u) ■ Vitdz — 4>i(u)b(x,u)dz + cj)i(u)g(x,u)dW + —cj)' 1 (u)G 2 (x,u)dt 

and since 

V1 M >{ = S u= ^Vu in V'(M. n x R) a.s. 
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we test the above against & € C™( R N ) and deduce 

d(l«>{, 0102} = ~{a(x, 0V1 U> 5, 0i02)dz - (b{x, £)£„=£, 0 i 0 2 )dz + (#(x, £)^=£ > 0i02)dIT 
+ ^{G 2 (x,£)£«=£, 0i02)dt 

= ~(a{x, ^)Vl u>i , 0 i 0 2 )d 2 + (6(x, 0 ^!u >6 0 i 02 )ds - (s(x, £)<%l u >{, 0 i 0 2 )dIT 
+ ^(<%(G 2 (x,£)<%l„>j),0i02)dt 

and conclude that the kinetic formulation (2.3) with the kinetic measure to = 0 is valid in the sense 
of V'($L n x R) a.s. In general, the kinetic measure is not known in advance and becomes part of the 
solution that takes account of possible singularities of the solution u. In particular, it vanishes in the 
above computation because we assumed certain level of regularity of u. 

In order to justify the formulation (2.5), note that, due to [5, Theorem 4] (similarly to Proposition 
5.2), if 0 £ G^(R W x R) then the composition cj>(9 t ) is the unique strong solution to 

d 4>(0 t ) + Vcj>(0t) ■ adz — d^(j){9t) bdz = 0, 

0(0o) = 0- 


As a consequence, if (F, to) solves (2.3) in the sense of ^'(R^ x R) a.s. 2 then applying formally the 
Ito formula to the product (F(t),<j)(0t)) we deduce 


( 2 . 10 ) 


where 


(F(t),(f>(9t)) = {Fo,(f>) — [ t (X7F.aA(0 t ))d Z + f {d^F b, 0(0 t ))dz 
Jo Jo 

- j\dtFg, m))dW + i j\ds(G 2 dsF), 0(0 t ))d t - m(l [o , t) c> £ 0(0 t )) 

- f (F,V<t>{0 t )-a)dz + [ (F, 9^(00 b)dz, 

Jo Jo 

f (VF-a,<f>(9 t ))dz + [ (dtFb,<l>(0 t ))dz- [ (F,V0(0 t )-a>dz+ f (F,9^(9 t )b)dz 
Jo Jo Jo Jo 

= - f <V[F0(0 t )],a)d* + f (0*[F0(0 t )],&)dz 
Jo Jo 

= [ ([F0(flt)],divo- d{b)dz = 0 

Jo 


due to the fact that diva — d^b = 0. Thus, (2.10) is a stronger version of (2.5) that does not require 
time dependent test functions. 

To conclude, let us also shortly discuss the connection with the notion of kinetic solution from [12] 
in the case of z t = t. As no rough integrals then appear on the left hand side of (2.3), the flow 
transformation method used in Definition 2.2 was not needed in [12, Definition 2]. Consequently, 
their notion of kinetic solution relied on solving the corresponding version of (2.3) directly in the 
sense of distributions. According to the above reasoning, this is possible in the case of sufficiently 
smooth solution to (1.1). Since the setting of [12] differs from ours in several ways: spatial domain, 
x-dependence of the coefficient A 1 integrability assumptions on a kinetic solution as well as on a kinetic 
measure, we will not compare the two notions of kinetic solution further. 


2 For instance ( F,m ) = (1 ,,>£, 0) from the discussion above. 
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2.4. The main result. Our main result reads as follows. 

Theorem 2.8. Let uo G L 2 (fl; L 1 (R N ))nL 4 (fl; L 2 (R N )). Under the above assumptions, the following 
statements hold true: 

(i) There exists a unique kinetic solution to (1.1). In addition, it satisfies 

<c(E\\uo\\h x +E\\u 0 \\ 4 Ll ). 

(ii) Any generalized kinetic solution is actually a kinetic solution, that is, if F is a generalized 
kinetic solution to (1.1) with initial datum l uo >{ then there exists a kinetic solution u to (1.1) 
with initial datum uq such that F = 1„>£ for a.e. (t,x,t;). 

(iii) If ui, U 2 are kinetic solutions to (1.1) with initial data u\, o and u 2 fi, respectively, then for 
a.e. t G [0, T] 

E||(ui(t) - u 2 {t)) + \\ L i < E||(ui, 0 ^ «2,o) + | \l\- 

2.5. Application to conservation laws with stochastic flux and forcing. It is immediately 
possible to apply our theory to conservation laws of the form (1.2), where B is a stochastic process 
that is independent on W and that can be enhanced to a stochastic process B for which almost every 
realization is a geometric p-Holder rough path. For instance, one may consider another Brownian 
motion but also more general Gaussian or Markov processes (the reader is referred to [18, Part 3] for 
further examples). 

Let B be a Brownian motion defined on a stochastic basis (Cl,P, (J^t),P) and let B(w) be a real¬ 
ization of its (Stratonovich) lift. For u> fixed, we set 

z = (z\z 2 ) := (B 1 ^), B 2 (w)) 

and Theorem 2.8 yields the existence of a unique kinetic solution 

u(u) G L 1 ( n ; L°°(0, T; L 1 {R N ))) n L 2 (fi; L 2 (0, T; L 2 (R N ))) 

together with the corresponding .^-contraction property. Let us now set 

ft = Cl x Q, P = P®P, p t = p t <&p t , P = P®P 

and 

B(lo,uj) := B(u>), W(u>,uj) := W(u>). 

Both these processes are Brownian motions on (Cl, IP , (Pt)- P) and if B denotes the (Stratonovich) lift 
of B then 

B(a), u>) = B(a)). 

Moreover, standard results in rough path theory guarantee that the RDE solutions to (2.6) and (2.7) 
driven by z = B(w,u;) coincide P-a.s. with the corresponding SDE solutions and in particular they are 
(J^t)-adapted. It remains to verify that the Ito integral in the definition of kinetic solution (2.5), which 
is now constructed for a.e. u> as an Ito integral on (f2, IP, (P t ),F), can be extended to ($7, P, (Pt), P). 
But that follows directly from its construction: the corresponding Riemann sums converge in P for a.e. 
ui, therefore by Fubini’s theorem they also converge in P. 

As a consequence, (1.2) or more precisely its kinetic formulation (2.5) is solved in the natural 
(stochastic) sense and Theorem 2.8 applies pathwise in Co. 


( 2 . 11 ) 


E sup |Kf)|| 2 i +E 

0 <t<T 


hd * 
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3. Elements of rough path theory 


In this section, we recall some basic notions and results from rough path theory which are used 
throughout the paper. For the general exposition we refer the reader to [18], [34], [35]. 

Let G[ p l(R d ) denote the step-[p] nilpotent free group over R d and let us consider the following rough 
differential equation (RDE) 


(3.1) 


d y = V (y)dx, 
2/(0) = 2/o, 


where V = (Vi,..., Vd) is a family of sufficiently smooth vector fields on R e and x : [0, T] —>• G^ (R d ) 
is a geometric Holder p-rough path, namely, it is a path with values in G^(R d ) satisfying 


d _ Ho1;[0,T] 


sup 

0 <s<t<T 


| 




< oo. 


We denote by C°’p“ Hol ([0, T]; G^ (R d )) the space of all geometric Holder p-rough paths. 

Definition of a solution to problems of the form (3.1) is based on Davie’s lemma (see [18, Lemma 
10.7]) that gives uniform estimates for ODE solutions depending only on the rough path regularity 
(e.g. p-variation) of the canonical lift of a regular driving signal x : [0,T] —> R d . As a consequence, a 
careful limiting procedure yields a reasonable notion of solution to (3.1). To be more precise, we define 
a solution to (3.1) as follows. 

Definition 3.1. Let x £ G°’p‘ Hol ([0, T]; G[ p l(R d )) be a geometric Holder p-rough path and suppose 
that ( x n ) is a sequence of Lipschitz paths with the corresponding step-[p] lifts denoted by x”, i.e. 


<-t :=S lp] (x n ) t = 1 


dx" 


dx” <g> ■ ■ • <g> dxl 


/o <v<t 


/0<vi<”*<i;[p] <t 


'[p] 


such that 

x" —> x 

uniformly on [0, T] and sup ra ||x n || A_ H61; r 0 T j < oo. We say that y £ G([0,T];R e ), also denoted by 
Tfiy)(0,2/o! x ), is a solution to (3.1) provided it is a limit point (in uniform topology on [0,T]) of 

{tt(v)(0, 2/o; x n )\ n > 1} 


where 7T(y) (0, yo; x n ) denotes the solution to the ODE 

d y n = V{y n )dx n , 


2/ n (0) = 2/0- 


Under a sufficient regularity assumption upon the collection of vector fields V, there exists a unique 
solution to (3.1) which defines a flow of diffeomorphisms. 

Theorem 3.2. Let p > 1. Let x be a geometric Holder p-rough path and suppose that the vector fields 
V = (Vi,... , Vd) are Lip 7+fc for some 7 > p. Then the following holds true. 

(i) There exists a unique solution to (3.1), say 7T(v)(0, 2 / 0 ; x) £ C([0,T];R e ), and the map 

<t> ■ ( t,y) £ [0,T] xR e £ 7 r (y) ( 0 , 2 /;x) t £ R e 
is a flow of C k -diffeomorphisms. 

(ii) There exists C > 0 that depends on |U|Li P 7 + fc an d ll x || i_ho 1 ;[o,t] such that for every multiindex 
a with 1 < | a | < k the following estimates hold true 3 

Sup \d a (f)(y) | !-H61;[0,T] < C, SUp | d a (j) 1 (2/) | i-Hol- [0 t] < C. 
yeR e p ’ ’ yeR e p ’ 


3 Here, | • li_H6i-[o t] d eno ^ es the Holder seminorm of a mapping taking values in R e . 
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Proof. The proof of these results can be found in [18, Proposition 11.11] and in [9, Lemma 13]. □ 

4. Contraction property 

Let us start with the question of uniqueness. Our first result gives the existence of representatives 
of a generalized kinetic solution that possess certain left- and right- continuity properties. This builds 
the foundation for the key estimate of this section, Proposition 4.5, as it allows us to strengthen the 
sense in which (2.5) is satisfied, namely, we obtain a weak formulation that is only weak in x, £ (cf. 
Corollary 4.2) and therefore using time dependent test functions is no longer necessary. 

Proposition 4.1 (Left- and right-continuous representatives). Let F be a generalized kinetic solution 
to (1.1). Then F admits representatives F~ and F + which are a.s. left- and right-continuous, respec¬ 
tively, in the sense ofV(M. N x R). More precisely, for all t* £ [0,T) there exist kinetic functions F* ,+ 
on Q x R^ such that setting F + {t*) = F* ,+ yields F + = F for a.e. {uj,t,x,£f) and 

(F+(t*+ e),0(0 t . +E )> —► (F+(f *),#?*.)> £4-0 V0€C , c 1 ( R ff xR) a.s. 

Similarly, for all t* £ (O,! 1 ] there exist kinetic functions F*'~ on LI x R w such that setting F~(t*) = 
F*’~ yields F~ = F for a.e. (ui,t,x,£) and 

(F~(f- e),0(0 t ._ e )> —►<.F-(O,^(0f)> U0 W> 6 C\ (R^ x R) a.s. 
Furthermore, there exists a set of full measure A C (0,T) such that, for all t £ A, 

(4.1) (F-{t),<t>{et)) = (F + (t),Wt)) V^C^xR) a.s. 

Proof. Step 1: First, we show that F admits representatives F~ and F + that are left- and right- 
continuous, respectively, in the required sense. It follows the ideas of [12, Proposition 8] and [22, 
Proposition 3.1], nevertheless, as our mixed rough-stochastic setting introduces new difficulties we 
present the proof in full detail. 

As the space Cf (R ;V x R) (endowed with the topology of the uniform convergence on any compact 
set of functions and their first derivatives) is separable, let us fix a countable dense subset V\. Let 
(f> £V i and a £ Cf (]0, T)). Integration by parts and the stochastic version of Fubini’s theorem applied 
to (2.5) yield 

T 

g^{t)a{t)dt + (F 0 , <j>)a(Q) = (m, d^(0 t ))(a) P-a.s. 

where 

(4.2) g^t) = (.F(t ), m)) + J o {dtFg, <j>(0„))dW j (d e (G 2 ^F, <j>{9 s ))d, s. 

Hence dtg<$> is a (pathwise) Radon measure on [0,T] and by the Riesz representation theorem g$ £ 
BV([ 0, T]). Due to the properties of RF-functions, we obtain that admits left- and right-continuous 
representatives which coincide except for an at most countable set. Moreover, apart from the first one 
all terms in (4.2) are almost surely continuous in t. Hence, on a set of full measure, denoted by f2^, 
( F(t ), 4>{6t)) also admits left- and right-continuous representatives which coincide except for an at most 
countable set. Let them be denoted by (F, </>(^)) ± and set flo = Since V\ is countable, Do is 

also a set of full measure. Besides, for cj> £ V\, (f,w) ( F(t,uj),(j)(9t)) ± has left- and right-continuous 
trajectories in time, respectively, and are thus measurable with respect to (f, oj). For <f> £ C* (R^ x R), 
we define (F(t, w), (^(Ot))^ 1 on [0,T] x D 0 as the limit of ( F(t,u>),(j) n (6t)) ± for any sequence (</>„) in V i 
converging to f> in the topology of C}. (R' v x R). Then as a pointwise limit of measurable functions 
(F(-,-),<j)(dt)) ± is also measurable in (f,w). Moreover, due to the uniform convergence (j) n to (f> and 
boundedness of F, it has left- and right-continuous trajectories, respectively. Consequently, let w £ flo 
and let J\f u C (0,T) denote the corresponding countable set of times, where 

(F(f,w),<(>(0t)} - ^ (F(t,w),</>(6' t )) + for some </> £ V L . 
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It follows from the Fubini theorem that there exists a set of full measure A C (0, T) such that for all 
t £ Ait holds true that t Jf u for a.e. w £ H, i.e. 

(F(t), m)y = (Fit), </>(0 t )) + e v 1 a.s. 

thus 

(4.3) = (F(t),m)) + V0eC' c 1 (K Jv xK) a.s. 

Let us proceed with the construction of F + , the construction of F~ being analogous. In the sequel, 
we will frequently write expressions of the type F(t, 7r t ) as a shorthand for the composition of F(t, •) 
with 7r t (similarly for F + etc.). Therefore, (F(t, n t ), </>) is understood as 

[ F{t,-K t {x,£))<l>(x,t)dxd£ 

J RW+i 

and since the flow 7r is volume preserving 

(F(t,TT t ),4>) = [ F(t,x,£)(l>(9 t (x,€)) dxd£ = (F{t),^{e t )). 

It is now straightforward to define F + (t ), t £ [0, T), by 

(4.4) (F + (t),<K0 t )) ■■= (F(t),m)) + , </> £ Cl(M. N x R) a.s. 

and observe that F + (t) is right-continuous in the required sense. Note that F + (t) is a.s. a well-defined 
distribution since, due to the flow properties of 7r and 9 , the mapping (j) K > <j>(8t) is a bijection on 
C^R^ x R) with the inverse mapping (j) i—>• </>(7r t ). Moreover, it belongs to L°°(M. N x R) a.s. Indeed, 
for oj £ fi 0 and </> £ T>\, if t £ [0, T) is such that 

(F{t^),m)) + = (F(t,^)A(0t)) 

then it follows immediately from the fact that F(t) £ L°°{ R w x R) a.s. For a general t £ [0,T) we 
take a sequence t n ~\ t with the property above and use lower semicontinuity. Moreover, seen as a 
function F + : Q x [0,T] —>• L^ oc (R N x R), for some p £ [1, oo), it is weakly measurable and therefore 
measurable. Hence, according to the Fubini theorem F + , as a function of four variables w,f, x,£, is 
measurable. 

Next, we show that F + is a representative (in time) of F, i.e. for a.e. t* £ [0, T) it holds that 
F + (t*) = F(t*) where the equality is understood a.e. in (w,x,£). Indeed, due to the Lebesgue 
differentiation theorem, 

1 r t *+ e 

lim — / F(t,K t (x,£)) dt = F(t*,Kr(x,£)) a.e. (u,t*,x,£) 

£-S>0 £ 

hence by the dominated convergence theorem 

lim — f (F(t, 7r t ), 4>) dt = (F(t*, 7r t .), </>) £ C^R^ x R) a.e. ( u>,t*). 

e->0 £ J t » 

Since the left hand side is equal to (F + (t*, 7r**), (f) for all t* £ [0,T] and u; £ Ho, it follows that 
(F + (t*),m*)) = (F(t*),H0t*)) W> £ C^(R n x R) a.e. (w,t*) 

which implies 

F + (t*)=F{t*) in L°°(R n x R) a.e. {u,t*) 

and Fubini theorem together with the measurability of F + and F regarded as functions of four variables 
ui, t, x, £ yields 

F + {t*,x,£) = F(t*,x,£) a.e. (u,t*,x,g). 

Step 2: Second, we prove that for all t* £ [0,T) 

(4.5) (F+(t*± £),</>) -^(F+(t*),<f>) 


V^eCj( R w xR) a.s. 
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Towards this end, we verify 

(4.6) {F+{t*+e,TT t *),<t>) —►(F + (i*,7r t .),^> V^C^R^xR) a.s. 

and observe that since the flow 7 r is volume preserving, testing in (4.6) by <j)(j !>) € Cg(M JV x M) yields 
(4.5). In order to prove (4.6), we write 

|(F + (f* +£,7 Tt-),(j>) - (F + (t*,7T t .),(/))| 

< | (F + (t* +£,7T t . ),(/)) - (F + {t* +£,TTr + e),<j>)\ + \{F + (t* +£,TTf +s ),<t>) - (F + (t* , 7T t . ), <f) \ . 

The second term on the right hand side vanishes a.s. as £ —>• 0 according to Step 1 and regarding the 
first one, we have 

\(F+(t* +£,7T t .),0) - <F+(r +£, n t .+e),<i>)\ = \{F+(t* +£),m*) ~ ^f+ S ))| 

and we argue by the dominated convergence theorem: according to Theorem 3.2, if supp<(> C Br, 
where Br C x R is the ball of radius R centred at 0, then cfi(9t *) and (j>(9t*+ e ) remain compactly 
supported uniformly in £ and their support is included in Br, ct i/ p . Besides, 

(a, 0) - 0 ) —► 0 V(x, ()eM"xM 

and since F + takes values in [0,1], the claim follows. 

Step 3: Now, it only remains to show that F + (t*) is a kinetic function on X = LI x R^ for all 
t* S [0,T). Towards this end, we observe that for all t* € [0,T) 

F n (t*,x,£) := — f F(t,x,£)dt 
£n Jt* 

is a kinetic function on X = Ll x l w and by (2.9) the assumptions of Lemma 2.5 are fulfilled 4 . 
Accordingly, there exists a kinetic function F*' + and a subsequence (n* k ) such that, on the one hand, 

F n *(t*)^F*’ + in L°°(Ll x R n x M). 

Since, on the other hand, we have due to Step 2 that 

F n *(t*) — >F+(t*) in V'{R n x R) a.s., 

we deduce that F + {t*) = F*’ + for all t* £ [0,T), which completes the proof. 

Step 4 : The proof of existence of the left-continuous representative F~ on (0, T] can be carried out 
similarly and (4.1) then follows immediately from (4.3). □ 

Remark that according to Proposition 4.1, a stronger version of (2.5) holds true provided F is 
replaced by F + or F~. The precise result is presented in the following corollary. 

Corollary 4.2. Let F be a generalized kinetic solution to (1.1) and let F~ and F + be its left- and 
right-continuous representatives. Then the couples (F + ,m) and ( F~,m ), respectively, satisfy for all 
t € (0 ,T) and every <j> £ C' 4 (M JV x R), a.s., 

(F+(t),m)) = (Fo,</>) ~ f (d^Fg,cf{9 s ))dW + \ [ (d ? (G 2 ^F>(0 s ))d S - m(l [o , t] 3^(0 a )) 

Jo z Jo 

and 

(F-(t),<K0 t )) = (F 0 , 0) - A dtFg, <K0,))dW + \ f (dei&dtF, </»(0 s ))ds - m( l m dtf(e B )). 

Jo z Jo 

Furthermore, setting F~( 0) := Fq, then for every <j) G x R) and t* G [0,T) it holds true 

(4.7) (F+(t*) -F-(t*),m*)) = -m{l {t . } dz<K0 t .)) a.s. 

^Note that we may assume without loss of generality that the cr-field T is countably generated and hence, according 
to [6, Proposition 3.4.5], the space L 1 ^) is separable. 
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In particular, if A C (0, T) is t/ie set of full measure constructed in Proposition 4.1, then for all t* £ A, 
m does not have atom at t* a.s. 


Proof. Consider (2.5) with a test function of the form (s,x,£) <->■ </>(x,£)a E (s) where <j> £ C* 
and 

r i, s<t, 

a E (s) = < 1 — t < s < t + e, 

[o, t + e < s. 

Due to Proposition 4.1 we obtain convergence of the left hand side of (2.5) as £ —> 0: 
r T 1 I t+e 

/ {F(s),(t>{6 s ))d s u e (s)&s = — / (F(s),(j>(0 s ))ds—> ~(F + (t), <j>(0 t )). 

Jo £ Jt 


Then we have 


m(a s {s)d^{9 s )) = f 

J \c 


d^4>(9 s ) dm + f (l —- —-W<£(0 S ) dm 

[0,t) xl w xR J[t,t+e)x i w xl ' £ ' 


L 


[0,t] xl N xR 


d^4>(d s ) dm 


and the convergence of the remaining two terms on the right hand side of (2.5) is obvious because they 
are continuous in t. Therefore, we have justified the equation for (F + , m). 

Concerning (F _ , m), we apply a similar approach but use the function 


{ 1 , S<t-£, 

t - £ < S <t, 

0, t < s 

instead. 

Finally, the formula (4.7) follows from (2.5) by testing by (t,x,f) i-» <j>{x, f)a e (t) where 

a e {t) = - min {(t — t* + e) + , (t — t* — e) - } 

and sending £ —> 0. As a consequence, m has an atom at t* if and only if F~(t*) ^ F + (t*) hence, in 
view of (4.1) the proof is complete. □ 


Lemma 4.3. Let F be a generalized kinetic solution to (1.1) and let F + be its right-continuous 
representative. Assume that the initial condition Fo is at equilibrium: there exists uq £ L l (Il x M. N ) 
such that Fq = 1 „ 0 >£. Then F + (0) = Fo and in particular the corresponding kinetic measure m has 
no atom at 0 a.s., i.e. the restriction of m to {0} x R w x R vanishes a.s. 

Proof. Let mo be the restriction of m to {0} x R^ x R. Then we deduce from (4.7) at t* = 0 (recall 
that F~( 0) = Fo) that 

(4.8) (F+(O)-l Uo>6 0) = -m o (%/>) V^C^xR) a.s. 

Let Hr be a smooth truncation on R such that 0 < Hr < 1, H r ( 0 = 1 if |£| < R and H R (f) = 0 if 
|£| > 2 R, \8^Hr\ < 1. For <f £ C}. (R ;V ) we intend to pass to the limit R —> oo in 

(4.9) (F+(0) - 1 U0>5 , (j>H R ) = -m 0 (^H R ) 

Since m and consequently mo is a finite measure a.s., the right hand side converges to 0 a.s. as R —> oo 
due to the dominated convergence theorem, whereas for the left hand side, we write 

(F + (0) - 1 Uo> £,(/>Hr) = (F + (0) - Io^AHr) - (l Uo >£ “ 
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We make use of Remark 2.6 and there introduced shorthand notation Xu 0 = Xi uo >t’ which together 
with the dominated convergence theorem implies the convergence 

(luo>£ 

Consequently, we deduce that also the remaining term from (4.9), i.e. 

(F + (0) — l 0 >f, (fiHn) 

is converging since all the other terms converge. Therefore, necessarily, 

<F+(0) - 1 0 >iAH R ) —► (F+( 0) - lo> £ ,0). 

As a consequence, we obtain 

/-F + (0,£) - lo>£d £ = u 0 a.s. 

Jr 

and using this fact one can easily observe that 

p(0“ f l«o>C-^ + (0,C)dC= f (luo>c-lo>c)-(^ + (0,C)-lo>c)dC>0 a.s. 

Indeed, p(— oo) = p( oo) = 0 and p is increasing if £ € (—oo, uo) and decreasing if £ £ (u o, oo). However, 
it follows from (4.8) that p = —too and since mo is nonnegative, we deduce that mo = 0. □ 

The following result allows us to restart the evolution given by (1.1) at an arbitrary time t* £ (0, T). 

Lemma 4.4. Let F be a generalized kinetic solution to (1.1) on [0,T] with the initial datum To. Then 
for every t* £ [0,T), t i —> F(t* + t) is a generalized kinetic solution to (1.1) on [0, T — t*] with the 
initial datum F~(t*). 

Proof. Let a £ C^([0,T)), </> £ C^(R N x R) and test (2.5) by (t,x,f) <f{x,f)a e {t) where 

! a(t), t < t*, 

a(i)( l-^T"), t*<t<t*+e, 


0, 


t* + e < t. 


In the limit e — > 0 we infer 

/* (F(t),<f>(fi t ))dta(t) d t - (T+(f),^ t .))a(f) + (F 0 , 0) 

(4.10) Jo 

= f (dzFg,(j)(dt))a(t)dW - ^ f (d i (G 2 d i F), f>{e t ))a{t) dt + m(l [0 ,t*]a5 ? (/>(6» t ))- 


Subtracting (4.10) from (2.5) and using (4.7) leads to 
[ {F(t),<K0 t ))dta{t)dt+ (F~(t*),m*)Mt*) 


= £ {d^Fg,(/)(dt))a{t)dW ~\J (G 2 9 ? F), <^(6» t ))a(t) dt + TO(l [t . iT] a5 ? ()>(6 l t )). 

Finally we take (t,x,£) i-A cf>(Tr t *(x,f))a(t) as a test function and deduce that 

J ( F(t ), 4{8 t ., t ))dta{t) dt + (F~(t*), <j>)a(t*) 

= J t (dtFg,<l>(8r, t ))a(t)dW - (^(G 2 ^F),^ t ., i ))a(t)di + m(l [t . iT] aa^(9 i ., t )). 
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As a consequence, t K > F(t* + t) is a kinetic solution to (1.1) on [0,T — t*] with initial condition 
F~(t*). □ 


We proceed with a result that plays the key role in our proof of reduction of generalized kinetic 
solution to a kinetic one, Theorem 4.6, and the lA-contraction property, Corollary 4.7. To be more 
precise, the relevant estimate that we need to verify in case of two generalized kinetic solutions Fi. F 2 
with initial conditions F\,o and Fyo, respectively, is the following: 


[ Fi(t)(l - F 2 (t))dxdf < [ Fi j0 (l - F 2j0 )dj:df. 

J R w xR J t"xt 

Indeed, setting F := F\= F 2 then leads to the Reduction Theorem and using the identity 

/ l« 1 (t)>|(l - l« 2 (t)>$) d£ = (ui(t) - u 2 {t)) + 

JR 

gives the L 1 -contraction property. However, it is not possible to multiply the two equations directly, 
i.e. the equation for F\ and the one for 1 — F 2 as it is necessary to mollify first. This is taken care of 
in Proposition 4.5. In order to simplify the notation, we denote F = 1 — F. 

In the sequel, let ( gs ) be an approximation to the identity on R w x R and (kr) be a truncation on 
R w such that kr = 1 on Br , supp^R C B 2 r, 0 < hr < 1 and IV^rr) < Rr 1 . 


Proposition 4.5. Let Fi, F 2 be generalized kinetic solutions to (1.1). Let s, t £ [0,T), s <t, be such 
that neither of the kinetic measures m\ and m 2 has atom at s a.s. Then it holds true 

E J Kr(x)(F+ (t), Qs((x, £) - d s ,t(-,-)))z,cr(F 2 (t),Q5{{x,t) - 0 S)t (-, •))>!/,C da; d£ 

-E J k r (x)(F-(s),0s({x,O - {■,■))) Z,cr {F 2 (s),e s ((x,£) - {■,-))) yX dxdf 

(4-11) < C(t - s) 1/p (e [ [ dm 1 (r,z,o-)+E f [d m 2 (r,y,()) 

\ J J(s,t] J J 

+ C(t - s) 1/p (eJ^ J \C\^ y{0 dydr + Ej s J \a\ 2 dn^ z (a)dzdr^j 

+ CRS(t - s) 1+1/p . 


Proof. Step 1: Applying Corollary 4.2 and Lemma 4.4 to the generalized kinetic solutions Fi, F 2 and 
to the modifier gs in place of a test function, we deduce that for all (a;,£) £ M w x R 

( F i(i),Qs((x,€) -6a,t(-r)))z,a 

(4.12) = & ~ (•’ ■)))*,* ~ J s ( d ° F ^9, 0 s((x, 0 - 0 s>r {-, -))) z a dW 

-mi(l[ s , t] (-)^^((a:,0 ~ 9 s , r (-,-))) (' G 2 d a F 1 ,d a g s ((x ,£) - 0 s>r (-, •)))„ ^dr 

and similarly 


( F t(t),Qs((x,£) -9 s ,t(-,-)))yX 

= {F 2 (s), gs((x, 0 - ('> -)))v ,C + [ ( d C p 2 g, g s ((x , £) - 0 s , r (-, ■))) c dW 

J S 

+ m 2 (l [sA (-)d ( g s ({x, f) - 0 s>r (-,-))) + ^ J (G 2 d c F 2 ,d ( gs({x,£) - 9 a>r (-,-))) yX dr 


(4.13) 
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J a(r)(G 2 d cr F 1 ,d (T gs{{x ,£) - 0 s>r (-, ■))) Z(T dr, 


[ a(r)(G 2 d ( F 2 ,d c g s ((x,£) - e Str (-,-))) y( dr, 

J S 


Our notation in the above is the following: the variable (:r,£) is fixed in both (4.12) and (4.13); in 
(4.12) we consider the test function (z,a) >-»• ps((x,£) — 9 s ,t{z,cr)) whereas in (4.13) we make use of 

(y,0 Qd((x,€) - Q s ,t(yX))- 

We denote 

Zi (t) ■= ~ f {d a F ig ,g s {{x,0 -9 8 , r {-,-))) dW, 

J S 

(t) := [ (d ( F 2 g,es((x,€) -^, r (v))) ( dWj 

J S ’ 

and observe that 
Mi : Cb([ s ) T]) -a M, 

a l-A 7711 («(•) d a Qs{{X, £) - 0 s , r (•, ’) 

and 

M2 : Cb([s,F]) — > R, 

a i-A m 2 (a(-) d c g s ((x,^) - 0 s , r (-, ■ 
are Radon measures hence 

t i-A Ml([s,t]), fe>/i 2 ([s,f]) 

are cadlag functions of bounded variation. 

Multiplying (4.12) and (4.13), we have 

E(F+(t), g s ((x , 0 - 0 s , t (-, •))>,,„(F 2 + (t), gs((x, £) - -)))w,C 

- W (S), 0«((x, 0 - (', -)))z,a(F 2 (s), 0 5 ((X, 0 - (', -)))y,f 
= E(^r(s),M5((a:,0 - (•, •))>*,^(i) + E(Ff (s), pg((x,£) - 0s(v)))z,<tM2([s.*]) 

+ EIi(t)(F 2 (s),p 5 ((x,C) - (•, -)))t/,C + EXi(t)X 2 (i) + IEXi(i)/^ 2 ([s, t]) 

+ EMi([s,t])(F 2 (s),q s ((x,£) - 6 s {-,-))} yX + Eg 1 {[s,t])l 2 (t) + E f i 1 ([s,t})g, 2 ([s,t}) 

= J\ + • ■ ■ + J&- 

For J5, we apply the Ito formula to the product of a continuous martingale and a cadlag function of 
bounded variation, that is I\(t) and g 2 ([s,f]), to deduce the following integration by parts formula 
(see e.g. [39, Chapter II, Theorem 33]) 

li{t)n 2 {[s,t]) = f M 2 ([s,r))dXi(r) + f I^r) d(i 2 {r) 

J S J (s,tl 


(4.14) 


which implies 


and similarly we obtain 


Next, 


J 5 = E f I\ (r ) dM2 (r ) 

Jr = e[ I 2 (r)dMi(r). 

Ji=E f (Fi(s),es((x,Q - {-,■))) z,a(dcF 2 g, gs({x,£) - 9 8>r (-,-))) d.W r = 0 

J S ’ 
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and similarly J 3 = 0. J 4 can be rewritten as follows 
J 4 =E((Ji(t),J 2 (t))} 


= -E 


J d(z, cr) ■ g(y, C)es((x, £) - 9 s , r (z, a))g s ((x, £) - 9 Str {yC))du^ z (a)dty 2 y {C)dzdydr. 


Next, we apply the integration by parts formula for functions of bounded variation (see e.g. [40, 
Chapter 0, Proposition 4.5]) to J 8 and deduce 

J 8 = E^i({s})/x 2 ({s}) +E f Mi([s,r])d/i 2 (r)+E / /x 2 ([s,r)) dm(r) 

J (s,£l J fs,£l 


— J $0 + J&l + J& 2 - 

Since the kinetic measures mi and m 2 do not have atoms at s a.s. due to the assumptions, it follows 
that J $0 = 0. Regarding J 8 i, we make use of the formula for /ii([s,r]) given by (4.12). Namely 

J 8 i=E / (Ff(r),Q 5 ((x,Z) - d s ,r(-,-))) z a dy 2 (r) 

J(s,t] 

-E(Ff (s), 0 s((x,€) - 9 S (-,-))) u 2 ([s,t]) - E [ I 1 (r)d/r 2 (r) = J 811 - J 2 - J 5 - 

Going back to the product of (4.12) and (4.13) we see that J 2 and J 5 cancel and for Jgn we write 


<^811 — E 


(s,t] 


Ff(r)g s ({x,0 - 9s, r (z,a))d c gs((x,(;) - 9 Sir (y, C))dm 2 (r, y, C)dzdcr 


- J Fl ( r ) G2 (yX)Qs{{x,0 - 9 Sj r{z,a))d c g 5 ((x,t;) - 9 s , r (y, C))dv>^ y (C)dydzd<rdr 

Let us now continue with J 82 . We apply the formula for /i 2 ([s,r)) which can be obtained from the 
formula for F~ in Corollary 4.2 using Lemma 4.4, cf. (4.13). It yields 

Js 2 = E [ (F 2 (r), g s ((x, £) - 9 s<r {-, -)))y£ dm(r) 

-E (F 2 (s),q s ((x,£) - (■,-)))y,CMi([s,i]) — E / l 2 (r)dm(r) = J S2 1 -J&- J 7 , 
where Jq and J 7 cancel and for J 82 i we have 


J 821 = —E 




F 2 (r)gs((x, £) - 9 s , r (y, C))d a gs((x , £) - 9 s , r (z, a))drrn(r, z, a)dyd( 


-E 

2 


J J F 2 {r)G 2 (z, cr)gs((x, £) - 9 s>r (y, ())d a gs{{x, £) - 9 s>r (z, a))dvl z (a)dzdydC,dr. 


Step 2: Finally, we have all in hand to proceed with the proof. Integrating (4.14) with respect to 
x, £ we obtain 

E J k r (x)F+ (t)F 2 ( t)g s ((x , £) - 9 Sit (z, <j))g s ((x, £) - 9 s , t (y , C))dird£dydCd 2 :der 

-E [ kr(x)F^ {s)F 2 (s)gs((x,£) - (z,a))gs({x,£) - (y, C))da;d£dydCd~dcr = h 4-b h, 


where 
I x = E 

h = -E 


(s,t] J 

f 

(s,t] 


KR(x)Ff(r)g s ((x,£,) - 9 a>r (z,a))d c g s ((x,^) - 9 a>r (y,())dm 2 (r,y,(), 

/ k r (x)F 2 (r)gs{{x,£) - 9 s , r (y,())d <T gs((x,€) - 9 B , r (z,a))dmi(r,z,a) 
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Above as well as in the sequel, when we integrate with respect to all variables r, x, z, er, y , we 
only specify the kinetic measures mi, m 2 and the Young measures v 1 , v 2 , the Lebesgue measure being 
omitted. 

Since the technique of estimation of I\ and / 2 is similar, let us just focus on I\. It holds that 

h = -E f [ k r (x)qs{{x,Z) - O s , r (z,a))g s ((x,€) - 6 » s>r ( 2 /,C))du r 1 ’+(CT)dm 2 (r, 2 /,C) 

J(s,t] J 

+ E f ( k r {x)F^{ r)d a Q S {{x,i) - 9 s , r {z : cr))gs((x,£) - 9 s , r (y,£))dm 2 {r,y,0 
J (s,t] J 

+ E f [ K R (x)Ff(r)Q S {(x,£) - 0 Str (z,a))d ( ps((x,^) - 8 s , r (y,())dm 2 (r,y,() 

J(s,t] J 

= III + 1 12 + Il3> 

where u 1,+ is the Young measure corresponding to Tj + and the first equality holds true because the 
sum of the first two terms is equal to zero due to 5 



which holds true for a.e. (w, z) and every (r, x, £). Consequently, In < 0. We will show that I 12 + 1 13 
is small if t — s is small. Towards this end, we observe that 


(4.16) 


d a es((x, 0 - 0 s ,r(z, a)) = -V x ,zes((x, 0 - 0 St r(z, a)) ■ d„ 0 8 tr (z, a) 


and 


J K R (x)V x ^g s {{x, £) - 9 a>r (z, a))g s ((x, £) - 0 a<r {y, C))dxd£ 


(4.17) 



J ^x^k r (x)qs{(x,^) - 9 S} r{z,a))gs({x ,£) - 0 a , r (y, C))dxd^ 


it follows that (recall that 9 X was defined in (2.7)) 



• ^x,iQs({x^) - 9 Sjr (z,a))gs((x,€) - 0 atr (y,())dm 2 (r,y,() 


(4.18) 



X Qs((x, 0 - 9 a ,r{z, a))gs((x, £) - Q a , r {y, C))dm 2 (r, y, £) 


1121 + 1 131 


J By (■, ■)fj we denote the duality between the space of distributions on R CT and ('}. (1^). 
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To proceed, we use the ideas from [20]. Namely, we observe that for every r £ [s, t], 9|O,(0 is Lipschitz 
continuous and its Lipschitz constant can be estimated by C(t — s) 1 ^ . Indeed, according to Theorem 
3.2 

sup |D 2 0 S ,O,01 = sup |D 2 (9 S ,0,0 — D 2 0 s ,sO,01 < C(r — s) 1//p < C(t — s) 1 ^. 

x,feR w xR x,£eR w xR 

Besides, we only have to consider the case of 

10,0 ^ 0,0, cr) < S, |(a;,0 - O,0, C)l < <5 

hence 

10,0,0 - O,0, 01 < 2<5 

and consequently 

(4.19) 10, cr) - 0,01 = I7r s , r (O,0> CT )) - 7T s , r (O,0, 0) I < C\9 s ,r(z : a) - 0,0, 01 < CS 

which implies 

(4.20) 


|0-0,0, cr) - 00,0,0 | < C(t - s) 1/p S. 


Let us now consider the remaining terms in the integral with respect to (x, 0 z, a) and use the fact 
that O, is volume preserving and <5|V^paK-) < Cg 2 s{-)- It yields 

J KflO) | Vx,105(0,0 - 0,0> CT )) | .5(0,0 - 0,0, OM^d-dcr 

(4.21) </|v„0 5 (0,0-0,0)1.00,0-0,0,0)d^d,d. 

C f c 

<j / .25(0,0 - 0, (7 )).00,0 - 0,0,0)dzd£dzd<7 = —. 

Plugging (4.20) and (4.21) into I'm we conclude that 

O 21 <C(t- s) 1/p E [ [ d m 2 (r, y, (). 

J(s,t] J 

Since due to Theorem 3.2 

( 4 . 22 ) sup o,oi = sup |O0*,O,O - O0 x , s 0>0l < C(t - S ) 1/p , 

(x,f)eR w xR (x,f)eR N xR 

we use the fact that |V x K/{(x)| < R~ 1 and that O, is volume preserving to deduce 

O 31 < CR~ l {t- s) 1/p E / [ 05(0,0 - 0,0, o')).a(0,0 - 0,0>C))dm 2 (r,2/,<) 

J 

< CR~ x {t - s) 1/p E f f dm 2 (r,y,() < C{t- s) 1/p E f f dm 2 (r,y,(). 

J (s,£] J J (s,£] J 

Let us proceed with I 4 and I 5 . Using the ideas from (4.16) and (4.17) we conclude 

O = J F 1 (r)G 2 {yX)nR{x)[d c 9 s , r {yX) ~ 00,0,0] 

' V x ,£05(0, 0 - 0,0, C)).5(0: 0 - 0,0, cr))dt' 2 y (0 
-^E^ J Fi(r)G 2 0 1 C)V I « fl 0)-U,Off) 


x . 5 ( 0 ,0 - 0 , 0 , 0 ) 05 ( 0,0 - 0 , 0 , 0 ) d ^, y (0 


+ 


Is I -O,O,o-))05(O,O - O,0,0)d0,(0d^, y (0 
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= /4I + I42 + A3 
and similarly 

h = f s J^2{r)G 2 (z,a)K R (x)[d a 9 s<r (z,a) -d c 9 Sir (y,Q] 

■ Vx,£Qs{{x,€) - O s> r{z,a))Bs((x,^) - 0 St r(y,C))dVr,z(, a ) 

+ ^ E / / F2 ( r ') G2 ( z ’ a ^ xKR ( x '> ' d ^ s A y ^) 

X Qs((x,0 - 0 St r(z,a))gs((x,£) - 9 s , r {y, C))du PZ (a) 

+ ^ E L I G2 ( Z ' a ^ KR ^ Q5 ^ X '^ ~ e s,r( z i a ))es(( x ,0 - 9s,r(y,())dVr,z( a ) d V?,y(() 

= hi + I 52 + ^ 53 - 

Using a similar approach as for /121 together with the fact that G 2 (yX) < C|C| 2 which follows from 
the assumptions on g , we deduce that 


hi + hi < C(t - s) 1/p (e j j G 2 {y , C)du r 2 y (C)dj/dr + J G 2 {z, a)d^ z (a)dzdr 
< C(t - s) 1/p (e J |CW, y (Odydr + Ejr j \a\ 2 dv\ z (a)dzdr 


and similarly to /131 we get 

I 42 +152 < CRT l {t - s) 1//p E J J \C,\ 2 dvl y {C,)dydr-\-CR~ l {t - s) l h^ J J \(j\ 2 diy pz (a)dzdr 
< C(t - sY^E J |(| 2 dv 2 y (C)d?/dr + C(t - s) 1 / p E J J |cr| 2 du^ z (cr)d 0 dr. 


Besides, 

h + I 43 + I 53 

r t 


= \ E j s J KR ( X ')\ 9 ( Z ’ <T ') ~ 9 (yX)\-Q s((x,£) - 9 St r{z,a))gs((x,£) - 9 a , r (y, C))du^ z (a)du 2 y (Q 

<C6 2 E f f e&{{x,h)-9sA z iv))Q5({x,C)-9sAyX))dvl' Z (a)dvly(C,) 

J s J Qr 

where Q R = B 2R . X xl^x B 2 R+ct i/ Pz xR„x B 2 R+CT i/ P y x and we used (4.19) and the Lipschitz 
continuity of g. Next, we estimate using (4.15), the fact that <5|Va;,^£?< 5 1(•) < Cg 2 s{-), (4.22) and the 
fact that 9 is volume preserving, as follows 

Qs((x,£) - 9 s , r (z, cr))du PZ (a)dz < / F^(r,z,a)d„g s ((x,^) - 9 S}r (z,a))dzda 

Vxl 

Fr (n 2, a)V x ,tQs{{x, £) - 9 s>r (z, a))d a 9 s , r {z , cr)dzdcr 


IB. 


2R+CT 1 /P ,z 


Jr n xR 


< C6 \t-s) l/p f Q 2 s((x,^) - 9 St r(z,a))dzda = CS |(t - s) 1//p 
JR N xR 


which holds true for a.e. (w,r) and every ( 2 ;, £). Hence 


h + h3 + h 3 < C6(t - s) 1 ^ 


s JB. 


2fi+CT 1 /P,y X ^ ^^2R,a:XR^ 


0s((x,O ~ 0 s ,r(y,C))daxl£du 2 y (C)d2/ 
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< CS(t - S 


cl*4(C)d V < CRS(t - s) 1+1 ' p 


5 2H + CTl/P,/ K C 

which completes the proof of (4.11). □ 

Finally, we have all in hand to prove the Reduction Theorem and the Z^-contraction property. 

Theorem 4.6 (Reduction Theorem). Let uq £ L 1 (~l L 2 (VL x R w ) and let F be a generalized kinetic 
solution to (1.1) with initial datum Fq = 1„ 0> £. Then it is a kinetic solution to (1-1), that is, there 
exists 

u £ L 2 (fl,L 2 (0,T-, L 2 (WL n ))) 
satisfying (2.4) such that F(t,x,f) = a.e. 

Proof. Let t £ [0,T), h > 0 and {0 = to < ti < ■■■ < t n = t} be a partition of [0 ,t\ such that 
| ti — ti +1 1 < h, i = 0,..., n — 1. According to Corollary 4.2, we can assume (without loss of generality) 
that the kinetic measure m has a.s. no atom at ti, i = 1 , ... ,n — 1 , that is, we only consider points 
from the dense set A C [0,T] which was constructed in Proposition 4.1. Besides, due to Lemma 4.3, 
m has no atom at 0 a.s. Furthermore, if v + denotes the Young measure corresponding to F + regarded 
as a kinetic function on Q x [0,T] x M^, then it follows from (2.9) that 


E / / / \£,\ 2 dv?,x(€) dxdt < oo. 

Jo J I N J R 

Consequently we may assume (without loss of generality) that for alli£{l,...,7i — 1} 


(4.23) 


E 


[ f l^ch'tUOd* < oo 
Jr n Jr 


and since zio £ L 2 (Ll x R w ) and u(f o x = S Uo ^ due to Lemma 4.3, the same holds true also for i = 0. 
Application of Proposition 4.5 yields 

(4.24) 

E J K R (x)(F + (ti ), g s ((x,£) - 9 ti _ ltU (■))) z ,a(F + {t i ), gs{{xA) - 9 ti _ uU (-, ■))) yX dxdf 

- E J kr(x)(F + (L- i), 0 S ((x, £) - -))) z , a {F + (U- 1 ), g s ((x, £) - (•, -)))y,C da;d £ 

< C/i 1/p E f f dm(r, y, C) + Ch 1/p E f ' [ |C| 2 dz/ r , y (C)dydr + CR5h 1+1 / p 

^ (ti — l ,ti] J J ti —1 J 

Now, we send 5 —>■ 0 first and then R —» oo. 

Regarding the first term on the left hand side of (4.24), since the flow is volume preserving we have 

{F + (ti),g s ((x,£) — ')))z,<r = (F + (t i ,Trt i _ 1 , u (-,-)),Qs((x,0 - (-, -))) z , a 

which converges to F + (ti, nt i _ lXi (x, £)) for a.e. (x,£) as <5 —> 0 and similarly for the term including 

F + (ti). Thus, if * = n we apply Fatou’s Lemma to estimate the liminf^o of the first term on the left 
hand side from below by 

E J K,R{x)F + (t n ,TTt n _ utn (xA))F + (t n ,TT tn _ u t n (x,f))dxdf 


= E 


J krW 


. (x, i))F + (t n ,X, £)F + (t n , x, 5)dxd^. 
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If i £ {0 ,... ,n — 1} we need to pass to the limit. To this end, we note that according to (4.23) and 
Remark 2.6 the functions Xf+(u) are integrable in and locally integrable in x. Therefore, we 
estimate 

(■ F + {U ), Qs((x,£) - 0 ti _ llti (-, ■)))z,*{F + (ti ), gs{{x,£) - 0 ti _ uti (-, -))) yX 
< f (F + (ti),Qs((x,0 - 0 ti _ uti {-,-))) z>(T , if f > 0, 

~~ \(F + (ti),Q S ((x^) -9 ti _ uU (-,-))) yi c, if ^ < 0, 

and observe that the left hand side converges a.e. in (w,x, £). Indeed, both functions F + (U) and 
F ( ti ) are locally integrable in (w,x,£) so their mollifications converge a.e. Next, we apply Remark 
2.6 to deduce that F + (ti) is integrable on x B 2 R+ct i/ Px x [0, oo) whereas F (ti) is integrable on 
Q x B 2 R+ct i/ p x x (—oo,0). As a consequence, 

(F + (ti), gs((x, £) ~ 0 ti _ lXi (-,-))) Z} a —> F + (ti,TT ti _ lA (x,£)) in L 1 (Xl x B 2 r x [0,oo)) 
and similarly 

{F + (U), es((x,£) -0 ti _! ,t 4 (v)))j/,C — >F + (ti,TT U _ 1 ,t i (x,€)) in x B 2 r x (-oo,0)) 

and for subsequences we obtain convergence a.e. in (ui,x,£). Therefore, we may apply the generaliza¬ 
tion of the Vitali convergence theorem [26, Corollaire 4.14] to deduce that 

(■ F+ (ti ), Qs((x,£) - -)))z,a(F + (ti), gs((x,0 - e ti _ uU (-, •)))«,c 

—> F + (ti,n u _ uti (x,£,))F + (ti,Tr ti _ ltti (x,£)) in L x (n x B 2 r x R) 

and accordingly to pass to the limit as 6 —> 0 in the first and second term on the left hand side of 
(4.14) to obtain terms of the form 

E J KR(x)F + (t i ,Trt i _ liU (x,£,))F + (ti,TT ti _ lM (x,$,))dxd £ 

= eJ KR(6 l_ lM (x,t))F + (ti,x,£,)F + (U,x,£,)dxdZ. 

In order to pass to the limit as R —> oo we intend to apply the dominated convergence theorem. To 
this end, it is necessary to justify that for all* = 0,..., n, 

F + (ti)F + (U) £ L 1 (I2 xR N x R), 

which can be proved inductively: if i = 0 then F + (0)F + (0) = 1 U0> ^(1 — lu 0 >c) = 0 hence the claim 
holds true. If it is true for some i — 1 where i £ (1,..., n} then the Fatou lemma together with the 
discussion above leads to 

E j F+ ( ti,x ’^ F+ ( ti,x, ^ dxd t ~ ,1™,® / KR ^U-i,t i ( x ^)) F+ ( t i-i^ x ^) F+ ( t i-i,x,^)dxd^ 

+ C/i 1/,p E f f dm(r, y, () + CTi^E f f \((\ 2 dv r ^y(C,)dydr 

J J J ti— l J 

= E J F + (t i - 1 ,x,£)F + (t i - 1 ,x,(,)dxd{, 

+ Ch l ^ p E f f dm(r, y, () + C/i 1/,p E / f \^\ 2 du rty (Qdydr. 

J (ti — l ,ti] J J ti—i J 

Thus Definition 2.1(h) and (2.9) give the claim for i. 

Altogether, sending 5 —t 0 and R —> oo gives for all* £ {1,..., n} 

E J F + (ti,x,£)F + (ti,x,£)dxd£ - E J F + (U-i, x, £)F + (fi_i, x, £)da:d£ 
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< Ch^PE / / dm(r, y, () + C7i 1/p E 


IC| 2 d^, y (C)d2/dr. 


Consequently, 

E J F+(t,x,Z)F + (t,x,0 dxd£-E J F+(x,OF+(x,(i)dxdt 
= E*/ F + {ti, x, £)F + (ti, x, £)dxd£ - E J F + (t i - 1 ,x,£,)F + (t i - 1 ,x,£)dxdf 

dm(r, y, C) + C'/i 1/p E J J |C| 2 d^ r ,y(C)dydr 


2=1 
< C7i 1/p E 


'(o,t] ■ 


and sending h —> 0 yields (using Definition 2.1 (ii) and (2.9)) 


E 


,£,£)-F + (i,£,£)dzd£ < E / F 0 (x,£)F 0 (x,£,)dxd£ = E / l„ o(x)>f (l - l„ 0 ( x)>f ) dxd^ = 0. 


Hence F + (t)( 1 — F + (t)) = 0 for a.e. (w,x,£). Now, the fact that F + (t) is a kinetic function for all 
t G [0,T) gives the conclusion: indeed, by Fubini Theorem, for any t G [0,T), there is a set E t of 
full measure in fl x such that, for all (uj,x) G E t , F + (t,x,£) G {0,1} for a.e. (gK. Recall that 
~d^F + (uj, t, x, •) is a probability measure on R hence, necessarily, there exists «:9x[0,T)x R^ —> R 
measurable such that F + (w, t, £, £) = l u ( w ,t >x )>£ for a.e. (uj,x,€) and all t G [0,T). Moreover, 
according to (2.9), it holds 

Eesssup / \u(t, x)\ da: = Eesssup / / |£| dvf x (£) da; < C 

0<t<T J rjv o <t<T Jrn J r 


and 



hence u is a kinetic solution. 


| u(t, a;)| 2 dxdf = E 




\£\ 2 dxdt < C 


□ 


Corollary 4.7 (.^-contraction property). Let u\ and U 2 be kinetic solutions to (1.1) with initial data 
ui,o and U 2 ,o, respectively. Then there exist representatives u\ and v ,2 of the class of equivalence ui 
and U 2 , respectively, such that for all t G [0, T) 

E|| (ui(i) - u 2 {t)) llj-i < E||(m 1j0 - u 2 ,o) + \\li■ 

Proof. First, we observe that the following identity holds true 

{ui-u 2 ) + = / l ui >£l« 2 >f d£- 

JR 

Let u\ and u 2 denote the representatives constructed in Theorem 4.6. Then proceeding similarly to 
Theorem 4.6, we apply Proposition 4.5 and obtain 

E|| (€ti (t) - w 2 (t)) + || L i =E J Ff(t,x,€)F2(t,x,g)d£dx 

— ® ^ ^-ui,o(x) >£-^-142,0 (a?) >£ -^H (^1,0 ^2,o) WlI 

which completes the proof. □ 
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5. Existence 


In the existence part of the proof of Theorem 2.8 we make use of the so-called Bhatnagar-Gross- 
Krook approximation (BGK), which is a standard tool in the deterministic setting. Its stochastic 
counterpart was established in [21] and could be understood as an alternative proof of existence to 
[ 121 - 

Let us now briefly describe this method. Towards the proof of existence of a kinetic solution to (1.1), 
which is (formally speaking) a distributional solution to the kinetic formulation (2.3), we consider the 
following BGK model 

d F e + VF £ -adz- d^F E bdz = gdW + j-d^G^F 5 ) d t + lu,>€ ~ F d t, 

( 0 . 1 ) AS 

F E ( 0) = F E , 

where the local density of particles is defined by 

u E (t,x)= f (F e (t, x,£) - l 0 >«)d£- 
Jr 


In other words, the unknown kinetic measure is replaced by a right hand side written in terms of 
F e . Heuristically, solving (5.1) is significantly easier and can be reduced to solving the homogeneous 
problem 

dX + VX ■ adz — d^X bdz = — d^XgdW + \-d^{G 2 d^X) df, 

( O.A ) A 

X{s) = X 0 , 


establishing the properties of its solution operator and employing Duhamel’s principle. The final idea 
is that, as the microscopic scale £ vanishes, F e converges to F which is a generalized kinetic solution 
to (1.1). The Reduction theorem, Theorem 4.6, then applies and as a consequence there exists u such 
that F = l u >£ and u is the unique kinetic solution to (1.1). 

We point out that the expected regularity in (x, £) of solutions to (5.4), (5.1) is low, namely L°°, 
and therefore the rigorous treatment of the above outlined technique requires a suitable notion of weak 
solution in the context of rough paths. Here the usual definition of distributional solutions encounters 
several obstacles due to rough regularity of the driving signals. Motivated by the discussion in Section 2, 
we introduce weak formulations which correspond to solving (5.4) and (5.1) in the sense of distributions 
but do not involve any rough path driven terms. 

In the sequel, we will use the following notation: for a € 1 we denote by Xa ■ ® t K the so-called 
equilibrium function which is defined as 


(5.3) 


Xa(£) — lo<£<c« la<f<0- 


5.1. Rough transport equation. This subsection is devoted to the study of the auxiliary equation 
(5.2). Our approach here as well as in Subsection 5.2 is rough-pathwise and therefore we work with 
one fixed realization u; and the probability space (fl,J^,P) remains hidden. The stochastic integral 
will then re-appear in Section 5.3 for the final passage to the limit. 

It can be seen that the Stratonovich form of (5.2) reads as follows 

dX + VX ■ adz — d^X bdz = —d^X g o dW + -d^Xd^G 2 df. 

Now we recall that (the stochastic process) A is the joint lift of z and W constructed in (2.1). Let us 
fix one realization A(w). This leads us to the study of the following rough transport equation 

( %X\ ( \d^G 2 b -g \ 

\XX ) \ 0 —a 0 ) 


dX = 


d(fA(w)), 
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However, for notational simplicity (and with a slight abuse of notation) we may rather write 

1 


(5.4) 


AX'- 


-ds-Xds-G 2 dt, 


XX ■ a dz — d^X b dz = —d^X g dw 

X{s) = X 0 , 

where w denotes the corresponding realization of the Stratonovich lift of W and the cross-iterated 
integrals between t, z and W are not explicitly mentioned. 

We introduce two notions of solution to (5.4). The first definition follows the usual approach from 
rough path theory which is based on the approximation of the driving signals (see [5, Definition 3]). 

Definition 5.1. Let (z, w) be a geometric Holder p-rough path and suppose that ( z n , w n ) is a sequence 
of Lipschitz paths such that 

S\p\(z n , w n ) = (z",w n ) —» (z,w) 

uniformly on [0,T] and 


sup | 


I — -Hoi < °°- 


Assume that for each n £ N 
dX n + VI" 


idz n - d e X n bdz n = -d(-X n gdw n + -d^d^G 2 dt 


t>N 


x n {s) = X 0 £ Cl( 

has a unique solution X n which belongs to C?([0, T\ x x R). Then any limit point (in the uniform 
topology) of ( X n ) is called a solution for the rough PDE, denoted formally by (5.4). 

Such a solution can be constructed by making use of the method of characteristics established in [5] 
provided the initial datum Xq is sufficiently regular. To be more precise, the associated characteristic 
system is given by 

^ dp? = —b(ipt) dz + g{ip t ) dw - i<%G 2 (p t ) dt, 

d tpf = a(ipt) dz. 

where p? and p® describe the evolution of the ^-coordinate and x-coordinate, respectively, of the 
characteristic curve. Let us denote by p Sjt (:r, £) the solution of (5.5) starting from (x,£) at time s. It 
follows from [18, Proposition 11.11] that under our assumptions p defines a flow of C 2 -diffeomorphisms 
and we denote by if) the corresponding inverse flow. 

Our first existence and uniqueness result for (5.4) is taken from [5, Theorem 4] and reads as follows. 

Proposition 5.2. Let Xo £ C?(R jV x R). Then there exists a unique solution to (5.4) given explicitly 
by 

X(t,x,£-,s) = X 0 (ip a>t (x,^)). 

We conclude that the solution operator S(t,s)X 0 = X 0 (ip s ,t(x,^)) is well defined on G?(R W x R). 
Nevertheless, as the right hand side makes sense even for more general initial conditions A'o which do 
not necessarily fulfill the assumptions of Proposition 5.2, the domain of definition of the operator S 
can be extended. In particular, since diffeomorphisms preserve sets of measure zero the above is well 
defined also if Xq is only defined almost everywhere. In this case, we define consistently 

S(t, s)X 0 = Xo^S'tiXjt;)), 0 < s < t < T. 

The aim is to show that this extension is the solution operator to (5.4) in a weak sense. Towards this 
end, it is necessary to weaken the assumption upon the initial condition and therefore a suitable notion 
of weak solution is required. 

In the following proposition we establish basic properties of the operator S. 

Proposition 5.3. Let S = («S(t, s),0 < s < t < T} be defined as above. Then 
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(i) for any p £ [l,oo], S is a family of linear operators on L P (WL N x R) which is uniformly 
bounded in the operator norm, i.e. there exists C > 0 such that for any Xo £ L P (R N x R), 
0 <s<t<T, 

(5.6) |[S(i,s)X 0 || L , <C\\X Q \\ L , 

x ’€ 

(ii) S verifies the semigroup law 

S(t , s ) = S(t, r) o S(r, s), 0<s<r<t<T , 

<S(s, s) = Id, 0 < s < T. 


Proof. The proof of (ii) as well as (i) in the case of p = oo follows immediately from the definition of 
the operator S and the flow property of if. In order to show (i) for p £ [1, oo), we observe that due to 
[18, Proposition 11.11] the map 

{s,t,x,£) i—!#«,*(£, 01, 

where J denotes the Jacobian, is bounded from above and below by a positive constant. Consequently, 
using the notation of Proposition 5.2, we have 

||*o||£, = [ jX(t,< / j Si t(z,<0;s)| , 'da;d£ 

x ’£ J l»xl 


K N xt 


I x(t, X , £ s)\ p I J if s ,t(x, 01 dx d£ > C\\X\\ p LP 


x,£ 


□ 


Definition 5.4. X : [s,T] x R w x K —>• R is called a weak solution to (5.4) provided 

dX(t, ifis.t) = 0 

( ‘ } X(s) = X 0 

holds true in the sense of T>'(lSt N x R), i.e. for all <f> £ C\( x R) and t £ [s, T] 

(x(t,ip a>t ),4>) = {x 0 ,4>). 

Corollary 5.5. Let Xq £ L° c (M. n x R). Then there exists a unique X £ L°°([s,T] x l w x R) that is 
a weak solution to (5.4). Moreover, it is represented by 

X(t) = S(t,s)X 0 . 


Proof. In order to prove existence, observe that if X is given by the above representation then clearly 
for all 0 < s < t < T 

X(t, (p St t(x,£);s) = Xo(ips,t°<Pa,t{x,£)) =X 0 (x,£) for a.e. (x,0 

and hence X solves (5.7). However, let us also give another proof of existence which in addition 
justifies Definition 5.4 as the appropriate notion of weak solution to (5.4). To this end, let us consider 
smooth approximations of Xq, namely, let (gj) be an approximation to the identity on x R and 
set Xq = Xq * gg. According to Proposition 5.2, there exists a unique X s which solves (5.4) with the 
initial condition Xq and is given by X s (t; s) = Xq (if s ,t) hence 

(5.8) X 6 (t,(p Stt (x,£)-,s) = X$(x,£). 

Moreover, 

X$^>X 0 in L°°( R w xR) 

SO 

X s (t;s) ^>X 0 (i/> s , t ) in L°°(R N x R) Vt £ [s,T] 
which justifies the passage to the limit in (5.8) and completes the proof. 
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Regarding uniqueness, due to linearity it is enough to show that any weak solution with Xq = 0 
vanishes identically. Let X be such a weak solution, i.e. it holds 

(X(t,'p a , t -,s),<l>)= 0 V^eC^R^xR). 

Testing by the modifier Qs((x,t;) — (•, •)) we deduce that for all (x, £) 

0= f X(t,<p att (y,C)\s) g s ((x,£) ~{y,C))dydC 
Jr n xR 

and since for a.e. (x, £) the right hand side converges to X(t, <p a ,t(x, £); s), the claim follows. □ 

Lemma 5.6. For all 0 < s < t < T it holds true that 

S(t, s)1q>j — 1q>{ = 0. 


Proof. It follows from (2.2) that for all x £ the solution to (5.5) starting from (x, 0) satisfies 
<p° s t (x, 0) = 0. Moreover, since the solution to (5.5) is unique, we deduce that 




> 0 , 

< 0 , 


if £ > 0, 
if £ < 0. 


Indeed, this can be proved by contradiction: let £ > 0 and assume that for some t £ [s, T) and x £ M. N , 
( Pst(. x >£) < 0. Since t i-A ip at (x,£) is continuous and <p® s (x,£) > 0 there exists r £ [s,t] such that 
ip° r (x,f) = 0. Now, we may take (p s>r (x,£) = (0, r (x, £)) as the initial condition for (5.5) at time r 
to obtain, on the one hand, 

<Pr,t(¥V (£,£)) = *(*>£) 

and, on the other hand, 

= FrA°’Fs,r(x,£,)) = 0 - 


which is a contradiction. As a consequence, S(t 7 s)1q>j = 1q>{ and the claim follows. 


□ 


5.2. Solution to the BGK model. Throughout this section we continue with our rough-pathwise 
analysis, that is, we consider one fixed realization of the driving path A(w) and therefore the underlying 
probability space (fl, J?,P) remains hidden. 

We will apply the auxiliary results for the rough transport equation (5.4) and establish existence 
and uniqueness for the BGK model (5.1). First of all, it is necessary to specify in which sense (5.1) 
is to be solved. As it can be seen in Definition 5.1, the solution given by Proposition 5.2 satisfies 
the equation (5.4) only on a formal level. This obstacle was overcome by the definition of weak 
solution (see Definition 5.4) which also permitted generalization to less regular initial data, namely 
Xq £ L°°(& n x R). We continue in this fashion and define the notion of weak solution to the BGK 
model similarly. 

Concerning the initial data for the BGK model (5.1), let us simply set Fq = l Uo >5 which is sufficient 
for our purposes. Nevertheless, our proof of existence and convergence of (5.1) would remain valid if 
we considered Fq = l«g>£ where («§) is a suitable approximation of uq. 


Definition 5.7. Let e > 0. Then F E £ L°°([0, T] x R^ x R) satisfying F e — lo>£ £ LHMxR^xR) 
is called a weak solution to the BGK model (5.1) with initial condition Fq provided 


dF £ (t,^ 0 ,t) 


° ‘fo ,« ~ <Po,t) df 


F e (0) = K 0 e 


(5.9) 


£ 
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holds true in the sense of V(U. N x R) 6 , i.e. for all (j) £ C^(M. N x R) and t £ [0, T] 

1 /"* 

(■ F e (t , (f 0 ,t), </>) = {Fq, <£) + -/ (l u e( 5)>c o y> 0 , fl - F e (s , ipo ,«), <^)ds. 

£ Jo 

The result reads as follows. 

Theorem 5.8. For any e > 0, there exists a unique weak solution of the BGK model (5.1) and it is 
represented by 

t 1 ft ± s 

(5.10) F e (t)=e~eS(t,0)F£ + - e ~S(t,s)l u e^ s)>i ds. 

£ Jo 

Proof. By Duhamel’s principle, the problem (5.9) admits an equivalent integral representation 

t 1 f^ts 

F e (t, <fio,t) = £-*Fq H— / e" 1 ^ £ l uE ( s ) >c ((po, s )ds 

£ Jo 

which can be rewritten as (5.10). Recall, that the local densities are defined as follows 

(5.11) u e (t,x)= [ (F e (t,x,£,) - 1 0 >{) d£ 

J R 

hence the function F e is not integrable with respect to £. For the purpose of the proof it is therefore 
more convenient to consider rather f E (t) = F e (t ) — lo>{ which is integrable and prove its existence. 
Moreover, we will show f E also admits an integral representation, similar to (5.10). Indeed, due to 
Lemma 5.6 and Corollary 5.5, 1 q>{ = S(t,s) lo>{ is the unique weak solution to (5.4) hence f e solves 


(5.12) 


, , e/ , v XV(t) ° Po,t ~ f e (t, ip o t ) 
d/ (t, <Po,t) = -—-d t 


f £ ( 0) = Xu 0 

in the sense of distributions. By a similar reasoning as above, (5.12) has the integral representation 

t 1 f^ts 

(5.13) f e {t) = e~ eS(t,0)xu o + - / e~ L ^ L S(t,s)xue(s) ds 

£ Jo 

and thus can be solved by a fixed point method. According to the identity 

[ |Xa - xp\ d£ = |a - /3\, a, /3 £ R, 

J R 

some space of ^-integrable functions would be well suited to deal with the nonlinearity term Xu s ■ Let 
us denote Jt? = L°°(0, T; L 1 (R JV x R)) and show that the mapping 

t 1 t 8 

=e~iS(t, 0)x„ o +- / e^ L ^S(t,s)xv( s ) ( is, 

£ Jo 

where the local density u(s) = f R ff(s,£) d£ is defined consistently with (5.11), is a contraction on Jif. 
Let g, gi, g 2 £ JFP with corresponding local densities v, Vi, v 2 . By Proposition 5.3 and the assumptions 
on initial data, we arrive at 

t 1 f^ts 

\\{JFg)(t)\\ L i < e“= ||5(t,0)xu o || L i +-/ ||«S(t, s)x v ( s ) || L i ds 




< C 


(l|Mo||ii+ sup \\v(s)\\ L i) 

V 0 <s<t J 


with a constant independent on t, hence 

< C'(||«o|| £i + ||ff||jr) < OO. 


®By ~\. u e (t)>^ 0 <po,t we denote the composition of (x,£) i —> l u e( t £C )>^ with (x,£) i—»■ 




30 


MARTINA HOFMANOVA 


Next, we estimate 


(JT g 1 ){t) - (JT g 2 ){t) 


1 rt t ^ 

1 — “ / ® £ ||<“*(^> 5 ) (Xvi (s) Xv 2 (s)) || 

x >* e J 0 

(j rt ^ 

— ~ I 6 £ ||Xvi(s) X”L»2 (s) 11 Z/ 1 

£ J o *•€ 

C ft t-s 

<— e _i ^||wi(s)-T 2 (s)|| L ids 

e JO 

C ft t s 

<— / e _i ^||gi(s)- p 2 (s)||li ds, 

£ Jo x ' 5 


ds 


so 

||— ^ 3 2 ||^ < C(l — e e )||gi — g 2 \\jtf 

and according to the Banach fixed point theorem, the mapping XL has a unique fixed point in Xf 
provided T was small enough. Nevertheless, since the choice of T does not depend on the initial 
condition, extension of this existence and uniqueness result to the whole interval [0,T] can be done by 
considering the equation on smaller intervals [0,T], [T, 2T\, etc. and repeating the above procedure. 
As a consequence, we obtain the existence of a unique weak solution to (5.1) that is given by (5.10) 
and the proof is complete. □ 


5.3. Convergence of the BGK model. In this final section we establish another weak formulation 
of the BGK model (5.1) that actually includes the stochastic integral and is therefore better suited 
for proving the convergence to (2.5). To be more precise, we prove the following result, which will 
complete the proof of Theorem 2.8. 

Theorem 5.9. Let f E = F e — lo>{- Then there exists u which is a kinetic solution to the conservation 
law (1.1) and, in addition, ( f e ) converges weak-star in L°°(Ll x [0, T] x x R) to the equilibrium 
function \u- 

We start with an auxiliary result. 

Proposition 5.10. F e £ L°°(Ct x [0,T] x R w x R) satisfies the following weak formulation of (5.1): 
let 0 £ C^(M. N x R), then it holds true a.s. 

(5.14) d<P E , 0(0 t )> = -(diF s g, 0(0 t )> dkF + i(d ? (G 2 S 5 F e ), 0(0 t )) di + i<l„.> £ - F E , 0(0 t )> d t. 
Moreover, it is progressively measurable. 

Proof. In order to verify (5.14), we test (5.10) by 0(0o,t) and obtain 

(5.15) (-F e (P0(0o,i)) = e _ =(5(t,0)Fop(6» o ,t)) + - [ e~ L ^(S(t, s)l u e (s)>? , 0(0 o , s (0 s ,t))) 

£ Jo 

According to Duhamel’s principle, it is enough to verify that if Xg £ L°°(M. N x R) then 

(X(t) :=S(i,a)* 0 ,#MM)) 

solves 

(5 16) d (X(t), <X0o , s (0 s ,t))) = -{diXg, 0(0o A9s,t))) dVF + i(S ? (G 2 d ? X), 0(0 o , s (0 s , t ))) dt, 

(X(s), 0(0 OjS )) = (X O ,0(0 O , S )). 

Indeed, let us calculate the stochastic differential of (F E (t), 0(0o,t)) given by the right hand side of 

(5.15) . Since according to the Ito formula applied to a product 

t 1 ft 

e -(X(t),0(0 o , s (0 s , t ))> = (X o ,0(0o lS )> - - / e ~(X(r), 0(0 o , s (0 s , r )))dr 

£ J ft 
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f 

J s 


l d (X(r),#MM)>> 


it follows due to (5.16) 

(F e (t),0(6> o ,t)) = (F 0 e ,^) - - [ e _ 5(5(r,O)F o e ,0(6»o,r-))dr 

£ Jo 

- f e-i (c>fyS(r, 0)F 0 %,0(0„,r)) dW r + ^ f e“i (^(G 2 9 ? [$(r, 0)F 0 e ]), 0(0 o ,r)> dr 

Jo ^ Jo 

+ -/ ( 1 K=( s )>{,</ , (^o. s ))-/ e _ ^(<S(r,s)l u6(s)>s ,(£(6fy s (6fy r )))dr 

£ Jo l £ Js 

- f e _z ^ £ (a 4 [5(r,s)l u£(s)>5 ]g,<(>(6io, s (6' S)7 .)))dW r . 

J S 

+ \J e ^ S )l« e (s)>{])) (t’{^o,s{0s,r))) dr ds 

= (F^cf>)- j\d^{r)g^^ r )) dW r + \j\d^G 2 d^{r))A{0o,r)) dr 

- - [ (F e {s),(j)(9o,s))ds+- f (lu e (s)>c j < / > (^o,s))ds, 

£ Jo £ Jo 


where the last equality is a consequence of deterministic and stochastic Fubini’s theorem and (5.10). 
Hence (5.14) is satisfied. 

Let us now justify (5.16). It was shown in [14, Theorem 8] that the RDE solution to the characteristic 
system (5.5) corresponding to the joint lift A constructed in (2.1) can be obtained as limit of SDE 
solutions to 


(5.17) 


d^’° = 9M) ° dVF - \dtG 2 {tf) dt - %?) dz n , 
dtf’ x =a(tp?)dz n , 


where the corresponding lifts (z™) approximate z as in Dehnition 3.1. Let ijj n and 9 n , respectively, 
denote the inverse flows corresponding to (5.17) and 


(5.18) 


d^° = -b(^)dz n , 
d7r"’ x = afy”) dz n , 


respectively. Let Xq G Cl(M. N x R) and (f> G G 2 (R W x R). Then setting X n {t) = X 0 (%/j™ t ) yields a 
solution to the stochastic transport equation corresponding to the characteristic system (5.17), which 
starts at time s from X 0 . Besides, <K#o s(^"t)) yields a solution to the transport equation corresponding 
to the characteristic system (5.18), which starts at time s from </>(9q s ). Therefore, applying the Ito 
formula to their product and integrating with respect to (x,£), we observe that the integrals driven 
by z n cancel due to the fact that diva — d^b = 0 (cf. a similar discussion in Subsection 2.3) and we 
obtain 

d(x n (t),m >a (0Z t ))) = -<^A" 5 ,^ 0 " s (^ t )))diT + i<a ? (G 2 a ? x"),^(^ s (^ t )))dt, 

(X (3)^(9 o, s )) = (Xo,m, B ))- 


As mentioned above, due to [14, Theorem 8], ip™ t (x, £) —»• ip St t(x, £) uniformly in t G [s, T] in probability, 
where i[> is the RDE inverse flow to (5.5). Due to invariance of the Lip 7 -norm under translation, we 
deduce (for a subsequence) that i/j n —>• tp uniformly in t G [0, T] and (x, £) G x R a.s. Indeed, since 
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ip™ t (x,£) is the SDE solution to (5.17) starting from (x,£) at time s, (p™ t := ip™ t (x,£) — (x,£) solves 

d =g((x,Z) + Tt) °dW - i<9 ? G 2 ((x,£) + (^) d t - b((x,£) + <p?) d z n , 
d<^”’ x = a((x,£) + <p?) d z n , 

<p n (s) = ( 0,0), 

and similarly (p Sit := y> Sj t(x, £) — (x, £) solves 

d<^t = ~b((x, 0 + fit) dz + g((x, 0 + <£ t ) dw - i<%G 2 ((x, 0 + (p t ) d t, 

d <Pt = a({x, 0 + <&) dz, 

0(s) = (0,0). 


Since for a family of vector fields V = (Vi, . .., Vd ) on R e it holds true that 

uny+oiiLip, = nn-)iiLip,, 

[14, Theorem 8] yields the convergence (p n —> (p uniformly in t £ [s, T] in probability which is indepen¬ 
dent of (x, £). Therefore (p n — > (p uniformly in t £ [s,T] and (x, £) £ R^ x R in probability and as a 
consequence, along a subsequence, ip n — > (p uniformly in t £ [s,T] and (x,£) £ R w x R a.s. Finally, 
according to Theorem 3.2 


supIC.fOs.O 

t,X,£ 


sup |C.tO^tfoO) - 

sup |(x, 0 - ip a , t (<p(x, £) + <p" it (x, 0 - y>(x, 0)1 

t,x,£ 


sup |^ s ,t(^ s ,t(x,0) + 

t,x,£ 


< sup |D-0 s ,i| sup |< t (x,0 - ¥? g)t (x,OI 

t,x,£ 

< Gsup |<p" t (x,0 - </? g ,t(x,OI -t 0 a.s. 

t ,X,£ 


and the claim follows. 

Moreover, 0 n —» 9 uniformly in s, t £ [0,T], (x,£) £ R w x R (see [5, Theorem 4]) and the same 
holds true for their first and second order derivatives with respect to £. Therefore, we may pass to 
the limit, apply dominated convergence theorem (for both Lebesgue and stochastic integral, see [39, 
Theorem 32] for the latter one) and (5.16) follows. If Xo £ L°°(R W x R) then we consider its smooth 
approximation Xq , apply the previous result and pass to the limit. 

The progressive measurability follows immediately from the construction. □ 


As the next step we prove a stochastic version of Proposition 5.3(i) for p = 1. 

Lemma 5.11. The family S = {S(t, s), 0 < s < t < T} consists of bounded linear operators on 
L x R w x R). In particular, if Xo £ L 1 (f2 x R^ x R) then 

(5.19) sup E sup ||5(t,s)X 0 || L i <GE||X 0 || l i . 

0 <s<T s<t<T x,t x,t 


Proof. Assume in addition that Xo is nonnegative, bounded and compactly supported in (x,£). In 
general, (5.16) holds true for all <j> £ G^R^ x R), nevertheless, since X £ L 1 (R Ar x R) a.s., which 
follows from Proposition 5.3, the assumption on the test function </> can be relaxed and we may take 
<f> = 1. Taking expectation, the stochastic integral vanishes due to the additional assumption and we 
obtain 

(5.20) E||$(*,«)* 0 || iiit <£11X011^. 
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Besides, the nonnegativity assumption can be immediately omitted by splitting Xq into negative and 
positive part. In the general case of Xq £ L 1 {SI x I w x R), we approximate Xq in L 1 (£d x R w x R) by 
X$ bounded and compactly supported such that 

E||X*||^<E||XoIL^- 

Apply (5.20) to Xq. Due to linearity of S(t, s) it implies that S(t, s)Xq is Cauchy in L 1 (D x R w x R), 
the limit is necessarily S(t, s)Xq so Fatou’s lemma yields (5.20) and the first claim follows. 

To obtain, (5.19), we test (5.16) again by </> = 1, take supremum and expectation. Applying 
Burkholder-Davis-Gundy’s and weighted Young’s inequalities and (5.20) we obtain 


E sup ||S(t,s)X 0 || L i < EHXolIzi +E sup 


s<t<T 


s<t<T 


(X, %/) 6W 
1/2 


<E||X 0 || iis +CE^ ||X||^ s dt 


<E||X 0 || l i +CE sup \\X\\ L i 




s<t<T 


1/2 


\ x \\ L ^dt 


1/2 


<E||Xo|| l i +Je sup ||X|| L 1 +ce/ T ||x|| l i dt 

1 s<t<T ’ 5 J s ,s 


<CE\\Xo\\l' +o E sup ||5(t,s)X 0 || L i 


2 


s<t<T 


and the proof is complete. 


□ 


Corollary 5.12. The family S = {S(t, s), 0 < s < t < T} consists of bounded linear operators on 
L 2 (D;L 1 (R JV x R)). In particular, if Xq £ L 2 (Cl\ L 1 (WL N x R)) then 

(5.21) sup E sup ||<S(t,s)X 0 || 2 1 <CE\\X 0 \\ 2 Ll . 

0<s<T s<t<T *■« 


Proof. The proof follows the lines of Lemma 5.11. The key observation is that if Xq is nonnegative, 
bounded and compactly supported in (&,£) then, on the one hand, 


E||«S(£, s)X 0 || 2 i <CE||X 0 || 


•CE 


[ (X&g) 

J S 


dW 


<CE||X 0 || 2 +CE /Vwilii dr 


and the Gronwall lemma implies 
(5.22) 


E||S(t,a)X 0 ||£i <GE||X 0 |||i 


On the other hand, Burkholder-Davis-Gundy’s inequality and (5.22) implies 


E sup ||<S(t, s)X 0 \f L i ^CEUXoll 2 ! +CE sup 

s<t<T x '* x s<t<T 


j\ x Ad) 


dW 


<CE||X 0 || 2 i 5 +CE 

<CE||X 0 |||i . 


\ x \\ 2 L idt 

x,£ 


□ 
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Proof of Theorem 5.9. In view of Proposition 5.10 we remark that F e also satisfies the following weak 
formulation of (5.9): let <f> £ C 2 ( x R) and a £ £?*([(), T)) then 

r T 


( F e (t ), <j>(e t ))d t a(i) d t + (F 0 e , </>)a(0) 


(5.23) 


= / (d^F E {t)g,cf(0 t )}a(t)dW(t) 



{d^G 2 d^{t))A{9t))a(t) d t 


~ F E (t),<j)(0t))a(t)dt 

£ Jo 

and a similar expression holds true for / £ , namely, it satisfies the weak formulation of (5.12). Taking 
the limit on the left hand side of (5.23) is quite straightforward. Indeed, according to the representation 
formula (5.10) it holds that the set (F e ) is bounded in L°°(f2 x [0, T] x x R), more precisely, F e £ 
[0,1], £ £ (0,1). Therefore, by the Banach-Alaoglu theorem, there exists F £ L°°{fl x [0, T\ x R w x R) 
such that, up to subsequences, 


(5.24) 

As a consequence, 


F E 


F in L°°(n x [0,T] x x 


(■ F E (t),(/)(9 t ))dta(t ) d t 


(^(G 2 %F £ (f)),<^))a(i) dt 


(F(t)A(0 t ))d t a(t)dt in L°°(fi), 
f (%(G 2 d £ F(t)),<£(0 t ))a(t)df in L“(f2), 


According to the hypotheses on the initial data, 

(F 0 £ ,^)a(0) = (l uo> {, </>)a(0). 

We intend to prove a similar convergence result for the stochastic term as well. Since 

<F e (t),^(^ t )))a(t)^<F(t),^( 5 ^0)Hi) in L 2 (n x [0,T]) 

and the stochastic integral <P i-a / Q T <PdW regarded as bounded linear operator from L 2 (fl x [0,T]) to 
L 2 (£i) is weakly continuous it follows 

[ T \d^F E {t)g,cf{B t ))a{t)dW{t)^ 


(d^F{t)g,^{e t ))a{t)dW{t) in L 2 (TL). 


Furthermore, since we established convergence of all the terms in (5.23) except for the third one on 
the right hand side, multiplying (5.23) by e gives the convergence of this remaining term to 0, that is, 

(5.25) [ — F e (t), <j)(8 t )}a(t) dt » 0 in L 2 (n). 

Jo 

As the next step, we show that the same convergence holds true if we replace the test function <f>(9)a 
by a general function /3 £ L 1 ([0,T] x R iY x R). To this end, recall that linear combinations of tensor 
functions of the form <pa , where (f> £ C 2 (M. N x R), a £ C^QO, T)), are dense in L 1 ([0,T] x R w x R). 
Hence if /? £ L 1 ([0, T] xt^xR) then there exists ^”=1 with </>* £ C^I^xR) and ay £ C£([0, T)), 
i = 1 ,..., n, such that 


{3 ^ ' Ayr? 


<<5 


J t,x,£ 


so, due to the fact that — F e £ [—1,1], we obtain 


[ F e (t),P(t,0 t ))dt <6 + ^2 [ (l^(t)>| - F £ (t),<l>i{0t))ai(t)dt 

Jo i=1 Jo 
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thus 

(!«'(()>{ - F £ (t),/3(t,e t ))dt — >0 in L 2 (tt). 

Consequently, we may take fi(t,nt) instead of /? to finally deduce 

l u e>£ — F E 0 in L°°{n x [0, T] x R N x M) 
and, in particular, for all <f £ C^(M), 7 

(5.26) (a € l u .> c - 0 in L°°(I2 x [0,T] x R N ). 

In order to obtain the convergence in the remaining term of (5.23) and in view of the kinetic formulation 
(2.5), we need to show that the term — F e ) can be written as d^rrf where m e is a random 

nonnegative measure over [0,T] x x R. bounded uniformly in e. If we define 

(5.27) m E {fi) = - f (l u .> c - ^(0) dC = \ [ (x»« (0 - 0) dC, 

it is easy to check that m e > 0 a.s. since F e £ [0, 1], Indeed, m e {— oo) = m £ (oo) = 0 and m e (t,x, •) is 
increasing if £ £ (—oo ,u e (t,x)) and decreasing if £ £ (u e (t,x), oo). 

Let us proceed with a uniform estimate for (u £ ) and (m £ ). 

Proposition 5.13. The set of local densities (u E ) satisfies the uniform estimate 

E sup |KC£)||ii < C'E||u 0 ||£i. 

0 <t<T 



Proof. It follows from the definition of u £ in (5.11) and from (5.13) that 


l {t) = e * S(t,0)xu o d£ + 


S(t,s)xu*(s) d^ds. 


Let us now define the following auxiliary function 


H(s) = 


S(t,s)xu*(s) d£ 


Then 


H(t) < e <= U(0) + (1 — e e ) max H (s) 

0<s<t 


and we conclude that H{t) < H( 0), t £ [0,T]. In order to estimate H( 0), we apply Corollary 5.12 and 
obtain 


E sup ||u e (t)|||i < E sup ||5(t,0)xu o ||r 1 < C'E||u 0 |||i. 

0 <t<T x 0 <t<T 


□ 


Proposition 5.14. For all t* £ [0,T] it holds true that 


(5.28) 


E 


I [0,t*] xl N xl 


dTO e (t, X, £) 


2 

+ E (f(t*),0 2 < C(E|| Uo ||1 5 +E||u 0 ||1i). 


^Here (•, -)^ denotes the duality between the space of distributions on R and Cc(R). 
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Proof. Here we follow the ideas of [20]. Let {0 = to < U < ■ ■ ■ < t n = t*} be a partition of [0,i*] 
with step size h, to be chosen below. It follows immediately from the formula (5.13) that for every 
ti £ [0,i*), t f e (ti + t) is a solution to (5.12) on [0,£* — U] with the initial condition and 

therefore the corresponding version of (5.14) holds true, namely, for all f £ C] (M iV x R) and all 
t £ [tii U-f-l] 

= -{dzf e g,(t>{Qt u t))&W + i(<9 ? (G 2 d 4 / £ ), <t>{6 ti ,t)) dt + (d ? m £ , <j>{6 tut )). 

Now, we need to test by </>(£) = £. Since f e £ L°°(0, T; L 1 (R Ar x R)) a.s., we can test by constants, 
in particular, we do not need compactly supported test functions. Therefore, let (j) R £ G X (R) be an 
approximation of <p which is bounded, monotone increasing, i.e. > 0, and preserves the sign, i.e. 

£sgn(/>#(£) > 0, and which satisfies |<9 |</>_r| < G uniformly in R. Using the weak formulation 

for f s above we deduce 

J (m £ ,5^fl(6» t °. iS ))d s + {f(t),<j) R (6^ t )) = () - J {d^f e g,cj) R (e 0 ti S ))dW s 

+ lj\d i (G 2 d^n,^el s ))ds. 

Observe that 0 < sgn(£)/ e (£) < 1 as a consequence of (5.13). Since also sgn(£)0° t (£) > 0 due to the 
assumption (2.2), the second term on the left hand side is nonnegative. Moreover, on a small time 
interval d^O®. t remains close to its initial value, that is, choosing h sufficiently small (which can be 
justified using Theorem 3.2) we may assume that for every i 


inf 

ti<t<ti -\-1 


d f &.. 


'Z 0 U,t ^ 


> - 


and consequently also the first term on the left hand side is nonnegative. Thus we deduce 


E 


(m £ , (d c 0 fl )((9° s )<%6»° s ) ds 


<CE(/ £ (b),fe) 2 + CE 


■e (nt),Mour 
2 


J (dtfg,M0°u, s ))dW s 


C E 


(dti&dtnM^ds 


Moreover, it follows from (5.13), Lemma 5.11 and Proposition 5.13 that 


E||/ e (t)||ii < sup E||<S(t, s)x u e (s) || 2 1 < GE||u 0 |||i 

0 <s<t<T x 


(5.29) 

and therefore the third term on the right hand side can be estimated as follows 


E 


(a ? (G 2 5 5 / £ ),^(0° jS ))d s 


<c{t-u) 2 n\uo\\ii- 


where the constant G does not depend on R due to the assumption on the derivatives of (j> R above. 
For the stochastic integral we have 

E f\dsrg,M0ls))dW s =Ef(3 { / £ j,fe(( s )) 2 d S 

ti Jti 

<CEj\f e , ds + G E £ ||r |||x e ds 

< GE [\ ft, M0° ti ,s)) 2 ds + C(t - u) ElKIlli 
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hence the Gronwall lemma yields 
(5.30) 

2 

(m e , (d(<f>R)(0f i>s )d^9f itS ) ds 


E 


E</ e (i), M0l, t )) 2 < Ch (e</ E (^), ) 2 + E||u 0 |||i). 


Therefore, if i = 0 then we estimate the first term on the right hand side of (5.30) by 


(5.31) 


E(f^^R ) 2 <E(xu 0 , 0 2 = ^n\no\\ 4 Ll 


and obtain by Fatou’s lemma 

rti 2 

(rrf, 1) dt ' 1B,/v2 


(5.32) 


E 


+ E(/ £ (t 1 ),6*o,ti) 2 < C h (E|| Uo ||i; +E||« 0 [| 


III 


In order to get a similar estimate on [ti,t 2 \ we go back to (5.30) and assume without loss of generality 
(using Theorem 3.2 again) that h was small enough so that for every i 

sup |£ — < 1. 

Consequently, by (5.29) and (5.32) 

E(/ E (ti),0 2 <CE(/ E (t 1 ) 1 ^) 2 + CE(/ E (t 1 ),(e-^J} 2 <C(E|| U o||i 1 +E|K||| i ). 

Iterating the above technique finitely many times, the claim follows. □ 


As a consequence of Proposition 5.13, the assumptions of Lemma 2.5 are satisfied for v^. x = S u er tx \^ 
and hence there exists a Young measure v tjX vanishing at infinity such that if —> v in the sense given 
by this Lemma. We deduce from (5.26) that d^F = —v hence F is a kinetic function. 

Next, we verify the second estimate from (2.9). Due to the definition of rrf in (5.27), it follows from 
(5.28) that 



e 



f e (t),0 d t 


+ E(f(t),0 2 <C{E\\u 0 \\l 2x +E\\u 0 \\ 2 Llx ) 


which implies 


(5.33) 



2 

<C(E|K||l2 x +E|M| i ). 


Rewriting the left hand side by the same argument as in (5.31) we deduce 


E 


K(t)||| 2 df 


< C^EIKIII; +EIML) 


and as a consequence the estimate (2.11) follows. 

Finally, in order to show that F is a generalized kinetic solution to (1.1), wc will prove that there 
exists a kinetic measure m such that, for all £ C'f(M JV x R) and a £ C^([0,T)), 

(5.34) f cj)(6t))a(t) dt -^4- f (d^m, (j)(9 t ))a(t) dt in L 1 (D). 

Jo Jo 

According to (5.28) and the fact that rrf > 0 a.s. we deduce that each rrf is a nonnegative finite 
measure over [0,T] x R w x R a.s. Besides, due to the convergence in (5.23), we obtain that for 
all (f> £ x R) and all a £ C,J([0,T)) the left hand side of (5.34) indeed converges weakly 

in L 1 (D) to some limit. Besides, due to Proposition 5.14, the set of measures (rrf) is bounded in 
I/)),(f2; A4&([0, T] x R w x R)), i.e. the space of weak-star measurable mappings from 12 to A4h([0, T] x 
R^ x R) with finite L 2 (f2; A4/,([0, T] x R w x R))-norm. Since L^(f2; A4&Q0, T\ x R w x R)) is the dual 
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of the separable space L 2 (D; Cb([0, T] x R w x R)), the Banach-Alaoglu theorem applies and yields 
existence of m £ L 2 (f2; A4{,([0, T] x R w x R)) such that, up to a subsequence, 

(5.35) m E —>■ m in L 2 (D; A4 b ([0, T] x R N x R)). 

This in turn verifies the convergence (and identification of the limit) in (5.34). It remains to prove that 
m is indeed a kinetic measure. Clearly, since all m E are nonnegative, the same remains valid for m. 
The points (i) and (ii) from Definition 2.1 follow directly from the construction of m and the uniform 
estimate (5.28). The remaining point Definition 2.1 (iii) can be justified as follows. Let (f> £ Co(R N x R) 
and define 

x E (t):=[ 0{x,t)dm E (s,x,£). 

xR 

If i? £ and 7 £ Z/°°(0,T) then by Fubini’s theorem 


E 

d [ 7 ( t)x £ (t)dt 

= E 

1 ? [ cj){x,Z)r(s)dm E (s,x,£,) 


. Jo 


. J[0,T]xR lv xR 


where T(s) = f^j(t)dt is continuous and -T(T) = 0. Since the right hand side converges to 


E 


' [0,T] xR w xE 


</>(x,{)r(s)dm(s,x,£) 


due to (5.35), we may apply Fubini’s theorem again to deduce that 


x(t) := / <f>(x,{)dm(s,x,£) 

I[0,t]xI N xR 

is a weak limit in L 1 (D x [ 0 , T ]) of progressively measurable processes and is therefore also progressively 
measurable. 

Altogether, we have proved that m is a kinetic measure and F is a generalized kinetic solution to 
(1.1). Since any generalized kinetic solution is actually a kinetic one, due to Reduction Theorem 4.6, 
it follows that F = l u >£ and v = S u , where u £ L 4 (fl; L 2 (0,T; L 2 (R N ))) is the unique kinetic solution 
to ( 1 . 1 ). 

The weak-star convergence of f E to Xu follows immediately from (5.24) and therefore the proof is 
complete. □ 
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